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ABSTRACT 


A method using electrical resistivity measurements to deter- 
mine the in situ porosity of marine sediments was investigated in 
the laboratory, and equipment for this purpose was designed, fabri- 
cated and tested, Formation Factor-porosity relationships deter- 
mined in the laboratory for three clays (kaolinite, illite, and 
montmorillonite), Providence silt, four sands, and four marine sedi- 
ments showed that porosity was predicted within + 2 percent. ‘The 
Formation Factors ranged from 1.1 to 5.9 while porosity ranged from 
26 to 93 percent. The particle size and distribution influenced the 
electrical resistivity of these sediments independent of porosity 
while particle shape did not. The laboratory equipment was econo- 
mical, safe, easy to operate and could be used to determine perme- 
ability, tortuosity and void ratio-log consolidation pressure of 
sediments. 

The Electrical Resistivity Measuring System designed for the 
Deep Ocean Sediment Probe consists of three interchangeable electrode 
arrays, the electronic circuit and the FM telemetry data link with 
6000 feet of coaxial cable. The predicted porosity values obtained 
with the System in the laboratory were within + 2 percent of the line 
of best fit obtained with the earlier laboratory equipment. Using 
the inner corer ring electrode oa the System was tested in the 


shallow water of Narragansett Bay and most of the data agreed well 





with the value obtained later with the laboratory equipment on the 
same sediment, 

The inner corer ring electrode array read continuous conducti- 
vity through the water/sediment column thus permitting a value of 

\ 

salinity to be obtained. The bottom water salinity obtained in this 
way was in good agreement with the interstitial water salinity of the 
cores. The system detected a core loss, touch down, and sample dis- 
turbance. This relatively low cost system allows accurate, quick 
in situ porosity determination and shows promise in other areas such 
as the monitoring of pollution, the prediction of other bottom sedi- 


ment properties and as a warning system on ocean structures, 
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PREFACE 


This dissertation is organized according to the standard plan, 
The text is divided into two sections. The first section includes a 
report of laboratory experiments on an indirect method of determining 
porosity of sediments and a review of the literature involving similar 
work, The instrumentation used in the laboratory is described in 
this section. 

The second section describes the design, fabrication, testíng, 
and evaluation of an in situ system used to determine porosity of 


marine sediments and the equipment used by others for similar work, 


xil 





THE USE OF ELECTRICAL RESISTIVITY TO PREDICT 


THE POROSITY OF SATURATED SEDIMENTS 


INTRODUCTION 


The purpose of this work is to investigate the accuracy of 
using an electrical resistivity method to predict the porosity of 
soil-water mixtures. Porosity in a two phase soil water mixture is 
the ratio between the volume of water to the total volume, The deter- 
mination of porosity is important to ocean engineering projects such 
as underwater acoustics, foundation engineering, submarine soil me- 
chanics and bottom survey work. At present, porosity is determined 
in the laboratory. If electrical resistivity is found to be an accu- 
rate method of predicting porosity, a practical method of in situ 
testing would be available. This method would be quicker, more eco- 
nomical, and allow more data to be taken than the present method, 

To verify the accuracy of the electrical resistivity method, 
extensive laboratory tests were conducted and the results were com- 
pared with theoretical and empirical work reported in the Literature. 
Most of the earlier work studying the relationship between electrical 
resistivity and porosity was conducted on formation rocks or disper- 
sions, The geophysical techniques used in well logging analysis 
involve electrical resistivity measurements in order to obtain: 

1. Subsurface correlation with formations in other wells, 

structural mapping. 

2. Lithology of the formations, 


3, Depth and thickness of productive zones, 





4, Formation porosities. 

5, Oil and gas saturation, 

6. Volume of oil in a reservoir, 

The porosity of formation rocks eee fromg3 to 20 percent, 
while the porosity of emulsions and dispersions ranges from 50 to 
90 percent. Although these ranges cover the porosity values of 
sediments (which range from 26 to 90 percent), extensive testing is 
necessary before it can be considered valid to apply the models and 
theories used for rocks, emulsions, and dispersions to sediments. 

Electrical resistivity has also been used in order to deter- 
mine the dielectric constant of pure crystals, the volume concentra- 
tion of certain biological materials such as red blood cells, and 
ground water detection, (Meredith?, Patten and Bennett?) 

Very little electrical resistivity work has been conducted on 
sediments, The limited work that has been conducted on sediments 
generally has not considered clay mineral content and type, particle 
size, shape, distribution and structure, and consolidation pressure, 

In this study the following parameters were varied in order to 
measure their influence on the prediction of porosity as determined 
by electrical resistivity measurements: 

1. Type of sediment (sands, silt, kaolinite, illite, and 

montmorillonite clays, and marine sediments), 

2, Size of particles, 

3. Shape of particles, 

4, Distribution of particle size, 


5. Structure of particles, 





6. Consolidation pressure, 

7. Salinity of interstitial water, 

The temperature was measured and a correction factor was applied in 
order to standardize all electrical resistivity measurements to 

25° Centigrade. By measuring the electrical resistivity as each 
parameter was varied, while at the same time determining porosity by 
the standard method, the influence of each parameter on the accuracy 
of this method of porosity prediction could be evaluated, 

Since the terms found in this work are used both in submarine 
soil mechanics and geophysics, they are defined below in order to 
eliminate confusion, 

The two-phase mixture of soil and water is considered a satu- 
rated sediment, or simply a sediment, The electrical resistivity 
measurements are taken on both the sediment and the interstitial 
water. The resistivity of the sediment is divided by the resistivity 
of the interstitial water to obtain a ratio called the Formation 
Factor. This ratio is used as a direct index of porosity. 

Electrical resistivity (specific resistance) is a measure of 
the ability of a solution, solid or mixture of unit length and cross 
section to resist the flow of electric current, In solids, current 
is carried by free electrons; іп solutions, the current is carried by 
ions. As electrons and ions move through solids and solutions, a 
frictional drag results. This frictional drag per unit volume is 
considered the electrical resistivity of the solid or solution, 


Electrical conductivity (specific conductance) is a measure of the 





ability of a solution, solid, or mixture, of unit length and cross 
section to carry an electrical current. Consequently resistivity 


and conductivity are inversely related measures of the same property. 


Studies dealing with formation rocks have related resistivity to 
porosity while studies dealing with emulsions and dispersions have 
usually related conductivity to porosity. 

Listed below are some of the symbols, assumptions, and defini- 
tions used in this paper: | 

dispersed phase (d) - solid soil particle, 

continuous phase (c) - interstitial water, 

mixture (m) - the mixture of the continuous and dispersed 

phases - sediment, 

The properties of these phases are assumed to be isotropic and uni- 


form, 


Em - resistivity of the continuous phase dispersed 
phase and mixture respectively. 

SOS - conductivity of the continuous phases dispersed 
phase and mixture respectively. 

In order to make relationships dimensionless the resistivity 


and the conductivity of the continuous phase is considered unity. 


These dimensionless notations are denoted by Roman capital letters: 


= 1, (1) 


YT 
d , 
Ка = , (2) 
C 





m 
Rn = т. , (3) 
С 
ko 
Ko D 146 (4) 
C 
k 
d 
Ka ig к. , (5) 
С 
Ка 
Е” (6) 
С 


Formation Factor (FF) which is often used in the Geophysical 
literature is equal to Rn and the reciprocal of к: 


1 
КЕ = Rn = 
m 


(7) 
Porosity is normally defined as the volume of the voids to the total 
volume, Porosity in a saturated sediment (two phase system) can be 
defined as the ratio of the volume of water to the total volume; 


Vo 
Porosity (п) = „IE, (8) 


total 


if saturated 


ker б 
= Т5 
total 





REVIEW OF THE LITERATURE 


The literature contains many theoretical and empirical models 
that Were developed in order to predict porosity from electrical 
measurement. The following review examines the models that might 
logically be applied to a saturated sediment. No one model can 
cover the wide range of sediment types and mixtures (sands, silts 
and clays). Factors such as particle size, shape, distribution and 
structure, and consolidation pressure influence model selection, 

In addition conduction of the clay fraction changes as the ionic 
concentration of the interstitial water varies. The first part of 
this review lists some of the properties peculiar to sediments that 
need to be considered in models in order to enable accurate porosity 
prediction from electrical measurement, 

Except for native metals, sulfides, a few oxides, graphite, 
and high-grade coal, the resistivity of most minerals that make up 
sediments is measured in millions of ohm-meters. Although the soil 
mineral is non-conductive, the soil particle carries an electrical 
charges The net electrical charge may arise from any one or a combi- 
nation of the following factors: 

1. Isomorphous substitution (important in illite and 

montmorillonite clays), 

2. Broken bond charges (important in kaolinite clays), 


3. Exposed hydroxyl groups, 


7 





4. Absence of cations in the crystal lattice, 

5. Presence of organic matter, 

The cation exchange capacity is a measure of the electrical 
charge carried by a soil particle. The magnitude of the cation ex- 
change capacity of a sediment depends on the amount and type of clay 
fraction present and the clay particle surface area, The clay frac- 
tion is usually defined as that part of the soil mass which has an 
equivalent Stokes diameter of less than two microns. Particles 
smaller than one micron are considered in the colloidal range. It 
is in the colloidal range that surface properties begin to dominate 
the physical and chemical properties of materials. 

For most clay particles the specific surface (the magnitude 
of the surface area per unit mass) is a good indicator of the rela- 
tive influence of electric forces on the behavior of the particle. 
Table 1 shows the cation exchange capacity and specific surfaces of 


the common clay minerals. 





TABLE 1 


THE CATION EXCHANGE CAPACITY AND SPECIFIC SURFACE 
OF KAOLINITE, ILLITE, AND MONTMORILLONITE CLAYS 


IR эол шшш сш ш. 


2 
meq. 
C100 Б) C) 


CLAY MINERAL —  —— — ты. 
CATION EXCHANGE CAPACITY” | SPECIFIC SURFACE 








Kaolinite 3 - 15 10 

1111ге 10 - 40 i 100 

Montmorillonite 80 - 150 1000 
E 


Jo : 4 
Lambe and Whitman 


When individual clay particles are dropped into water the ions on 
the surface of the particle tend to move out from the surface to be- 
come exchangeable ions. Water molecules are attracted to these ions 
and the two together form what is called the double layer around the 
clay particle. The thickness of the double layer is defined as the 
distance from the soil surface to the point where the ion concentra- 
tion is in equilibrium with the interstitial water, 

Many factors influence the size of the double layer, The 
Gouy-Chapman Double Layer theory is often used to explain the double 
layer and the factors affecting it. This theory assumes that there 
is no interaction between clay particles and that the clay particles 
are plates. This theory states that the size of the double layer is 


directly proportional to the temperature and dielectric constant of 
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the solution and inversely proportional to the concentration and 
ionic valence of the electrolyte in the interstitial water and the 


electric charge on the clay particle. Other factors influencing 


the double layer are рН апа size of ions in solution, 

If the cation exchange capacity and the specific surface are 
low and a small double layer exists, almost all the interstitial 
water would be free water and would not be attracted to the mineral 
skeleton, Conversely if the cation exchange capacity and the specif- 
ic surface are high and a large double layer exists, it would be 
quite difficult to separate the clay mineral and interstitial water 
phases. 

One of the earliest works examining the effect of charged 
clay particles and the associated double layer on electrical resis- 
tivity was done on reservoir rocks by Winsauer et = These in- 
vestigators concluded that charged clay particles do influence the 
conductivity of shaly sand and limestone formations, but that this 
influence became negligible at high ionic concentrations of inter- 
stitial water. 

After conducting laboratory experiments, Hill and Milburn® 
made similar conclusions and developed empirical equations for 
reservoir rocks relating Formation Factor to the resistivity of the 
interstitial water, porosity and the cation exchange capacity. In 
their equations as the electrical resistivity of the interstitial 
water decreased, i.e., became more saline, the effect of the cation 
exchange capacity decreased, 


One of the first works to study the effect of the cation ex- 
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change capacity on the electrical resistivity of a sediment was con- 
ducted by posse Using kaolinite clay and sodium chloride solutions, 
he showed that the Formation Factor remained constant when the water 
resistivity was less than 80 ohm-cm, However when water resistivity 
was greater than 80 ohm-cm, the Formation Factor decreascd, Water 
that has an electrical resistivity of 80 ohm-cm would have a salin- 
ity of seven parts per thousand, Sea water having a salinity of 

35 parts per thousand would have a resistivity of 25 ohm-cm, This 
work agrees with the work on reservoir P The effect of the 
charge of kaolinite clay particles on the electrical resistivity ap- 
peared to be negligible at high ionic concentrations of interstitial 
water, 

Sheeler et A expanded the knowledge on sediments by study- 
ing several types of clay. He found that for kaolinite clay the ef- 
fect of cation exchange capacity on the electrical conductivity and 
thus Formation Factor persisted up to a 3 percent sodium chloride 
solution, then became negligible. However when using montmorillonite 
the ion effect persisted up to an 8 percent sodium chloride solution, 
Sea water is approximately 3.5 percent sodium chloride. This study 
suggests that the type of clay present in a sediment such as 
montmorillonite may influence the electrical resistivity measurements 
due to the cation exchange capacity even though the salinity is 
higher than sea water, 

Boyce” reported that the percentage of clay size material 
varied inversely with the electrical resistivity of fine grained. 


marine sediments obtained from the Bering Sea. Boyce explained these 
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results could be due to: 


l. a specific conductivity increase resulting from a 
double ionic layer of clay mineral in the clay size 
fraction. 

2. a secondary conductivity increase caused by an increase 
in porosity, since the clay size fraction and porosity 
are directly related. 

Fritsch and Tauber s studied kaolinite, illite, and 
montmorillonite clays at low and high salt concentrations, They 
observed that the electrical resistivity measurements were over- 
whelmingly influenced by the specific surface area of the clay par- 
ticles when distilled water was used, but at high salt concentrations 
the resistivity measurements were independent of this influence. As 
illustrated in Table |, the specific surface area is an indicator 
of cation exchange capacity. These authors concluded that further 
work was necessary in order to establish the limits of this influ- 
ence for different clays at different concentrations, 

Atkins and smith studied sodium montmorillonite, calcium 
montmorillonite and kaolinite, which have cation exchange capacities 
of 74.0, 60.1 and 16.1 meq/100 g respectively. They found that the 
Formation Factor ‘remained constant even though the interstitial 
water conductivity ranges were l.l to 8,7, 0.1 to 19.7, and 0,009 
to 8,5 mho/meter respectively. They found the same result with 
illite clay, but the range of interstitial water conductivity used 
was not given, 


Kermabon et al, |4 further studied the effect of ionic solu- 
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tions on the cation exchange capacity of clay particles. In their 
laboratory experiments the marine clay particles had a grain size 
smaller than 3.9 microns and the sodium chloride solutions ranged 
from 0.7 to 10 percent, The porosity of the clay sediment samples 
was constant. Throughout this wide range of water salinity the re- 
lation between the conductivity of sediment and the conductivity of 
the interstitial water remained constant; therefore, any influence 
the charge of the clay particles had on conductivity was masked by 
the ionic solutions, 

ГЕН? found similar results in laboratory experiments with 
"coarse" clay. He used sodium chloride solutions ranging from 
0.5 to 10 percent. 

The results of these studies indicate that in fresh water the 
charge on clay particles (cation exchange capacity) greatly influ- 
ences the conductive properties of the sediment, This charge and 
consequently its effect on electrical resistivity measurements in 
fresh water depends on the amount and type of clay present. It is 
apparent that in a fresh water environment cation exchange capacity 
must be accounted for by the model or by empirically measuring and 
correcting for its effect on resistivity measurements if accurate 
porosity predictions are to be made. However in a marine environ- 
ment, the influence of the cation exchange capacity of clays on 
electrical resistivity measurements appears to be negligible in most 
cases. The reason for this difference between fresh and salt water 
can be deduced by visualizing the current flow in solution as the 


migration of free ions. In a sediment of clay and fresh water, most 
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of the ions in solution are under the influence of clay particles 
and therefore the current flow through the sediment is inhibited. 
If the interstitial water is saline, there are sufficient ions to 
neutralize the clay particles and allow the current to flow freely. 
If the clay's cation exchange capacity is high, more ions will be 
required to neutralize its effect on conductivity. Although cation 
exchange capacity is usually neglected in the marine environment, 
more testing is needed before this influence on electrical resisti- 
vity measurements for all types of clays can be assumed to be 
negligible. 

The above studies indicate that the cation exchange capacity 
of some clays need not be considered as an influence on electrical 
resistivity measurements in the salinity of ocean water. However 
other characteristics of clay must be considered if accurate poros- 
ity predictions are to be made from electrical resistivity measure- 
ments, The clay type, sand and silt present in a sediment determine 
the size, shape, and distribution of particles, and influence parti- 
cle structure, 

The M.I.T. soil classification of sediment particle size! 
is given below: 

from 0.06 mm to 2,0 mm is considered sand, 

0.002 mm to 0,06 mm is considered silt, and 

a particle smaller than 0,002 is considered clay sized. 

The particle distribution of a sediment is defined as a measure of 
spread such as standard deviation from the mean. If the results of 


past work are applied to sediments, it would appear that the parti- 
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cle size and distribution might influence electrical resistivity 
in the following ways: 

If the porosities, particle shapes, and standard deviations 
of two distributions are the same for samples with different mean 
particle sizes, the electrical resistivity will usually be lower 
for the larger particle size. This effect was suggested in experi- 
ments with glass beads by Meredith! and by Wyllie and . 2. 

As the particle size increased, the electrical resistivity de- 
creased, 

If the porosity aad nean particle SR and size are the same 
while the distributions of two sediments are different, the electri- 
cal resistivity will be higher in the'sediment with the greater as- 
sortment of particles. This effect was suggested by studies on rocks 
by Semenov"? and work conducted on glass beads by Meredith.. 

The effect of shapes in a sediment is complicated, A sedi- 
ment may be made up of particles of different shapes ranging from 
spheres to plates. Usually particles in the sand and silt size range 
are closer to a sphere in shape than the plate-like clay particles. 
The different particle shapes affect the path lengths that migrating 
ions take when current flows. If two sediments have approximately 
the sane particle size and spread and equal porosities but differ in 
particle shape, the electrical resistivity would be expected to be 
higher for a sediment composed of disc-shaped particles than a sedi- 
ment made up of spherical shaped particles. This effect was sug- 
gested by studies on various shapes by Wyllie and бс А Discs, 


cubes, cylinders, and triangular prisms resulted in higher electri- 
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cal resistivity values than spheres at the same porosity values. 
They also found that Ottawa sand (sphericity = 0.895) and glass 
beads (sphericity = 1.00) had similar Formation Factor versus poros- 
ity relationships. 

Meran, found that the conductivity of dispersions at a 
porosity of 50 percent resulted in a 10 to 20 percent error if the 
shape difference between spheres and discs was not taken into con- 
sideration. The error decreased with increasing porosity. 

Sana theoretically developed equations for formation 
rocks based on shape. Semenov's equations £o rocks of spherical 
shaped grains and disc shaped grains (having the same distribution) 
resulted in Formation Factor values of 2.7 and 10.1 respectively at 
a porosity of 50 percent. 

Atkins and ced stated that at the same porosity angular 
particles have higher Formation Factors than spherical shaped par- 
ticles. A shape factor was determined for five clays and one mica, 
Shape was inferred to be the major factor causing differences be- 
tween the clays. No electron microscope photographs were taken, and 
data such as mean particle size and distribution were not given. 

In this paper, particle structure will refer to the orienta- 
tion of individual soil particles in the sediment mass. Particle 
structure in sands and silts usually depends on packing. The type 
of packing is important because it affects the porosity and the path 
length required for migrating ions aad consequently the electrical 
resistivity measurement. The two extreme types of packing of uni- 


formly sized spherical particles are loose or cubic packing and 
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dense or rhombohedral packing as shown in Figure l, 


FIGURE 1 


SPHERICAL PARTICLES IN A DENSELY AND A 
| LOOSELY PACKED CONFIGURATION 


CUBIC PACKING RHOMBOHEDRAL PACKING 
(loose) (dense) 

Structure in clays is greatly с. by those factors af- 
fecting the size of the double layer. In clays the net electrical 
force between particles is the result of many repulsive and attrac- 
tive forces. If the double layer of clay particles is small, the at- 
tractive forces are dominant and the clay particles will move toward 
each other and become attached (Flocculant Structure). If the double 
layer is large, the repulsive force is more dominant and the clay 
particles will tend to move away from each other (Dispersed Struc- 


ture). The two clay structures are shown in Figure 2. 


FIGURE 2 


ORIENTATION OF CLAY PARTICLES IN A 
FLOCCULANT AND A DISPERSED STRUCTURE 
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FLOCCULANT DISPERSED 
STRUCTURE STRUCTURE 
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In the flocculant structure, large voids may be enclosed by 
the clay skeleton and result in sediments of very high porosities. 
On the other hand the dispersed structure results in lower porosity 
values, Aside from the facts of porosity differences, it would ap- 
pear by visualizing the two structures that the flow of current 
through these two structures would also be different, 

If the current path becomes more tortuous its length is 
increased,  Tortuosity of the current path may be defined as the 
ratio of the average statistical length of current flow lines to the 
simple physical length in a straight line from one end of a sediment 
sample to the other, 

By definition the resistance (R) between flat faces of a 


cylindrical sediment sample is: 


zr —© : (10) 


where: 
] is the length of the sediment sample, 
A is the cross-sectional area of the sediment sample, 
la is the effective average statistical length of the sedi- 
ment through the pore structure, 
A, is the effective average statistical cross-sectional area 
of the pore structure. 


In terms of Formation Factor: 
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where t is the tortuosity factor previously defined. Porosity can 


be expressed as follows: 


as qom As 
т d AA teo (12) 


If the porosity is the same equations (1) and (12 can be combined 


and tortuosity may be defined as a function of porosity and Formation 


Factor: 
к = ИРИ, (13) 
or 
t2 
FF = 8 (14) 


As the structure of a sediment becomes more complex and the 
current paths become more devious, the electrical resistivity will 
increase. It is at this point that the knowledge of both the elec- 
trical and hydraulic properties of sediments and the relationship be- 
tween them becomes important if physical properties are to be proper- 
ly interpreted by electrical resistivity measurements. ¡Much work has 
been done on the flow problem through porous media. The general flow 


formula, be the flow fluid, electrical, thermal or chemical, is: 


F = L*X+As . (15) 





20 


F 15 the flow rate, 

X is the gradient, 

A is the cross-sectional area, 

L is a coefficient of conduction, 

Numerous equations directly inter-relating hydraulic, thermal, and 
electrical properties of porous, solution-saturated materials have 
been derived, For example, Hutt and p developed prediction 
equations by investigating thermal and electrical conductivities of 
sandstone rock and ocean sediment and relating the two conductivi- 
ties through the common parameter of porosity. Although such a sys- 
tem seems reasonable it is important to remember that the success of 
any electrical method for measuring thermal conductivity in a sedi- 
ment would depend on the range of electrical conductivity of the 
soil particle encountered in a given situation. Therefore it is 
imperative when inter-relating flow problem that these differences 
be understood, 

Schopper contends that if quantitative evaluations between 
electrical and hydraulic properties are to be improved, the differ- 
ences between electrical and hydraulic tortuosity should be taken 
into consideration, The essential difference is that hydraulic 
tortuosity takes into account a direction factor in addition to the 
normal tortuous path length, while the normal tortuous path length 
is the only factor involved in electrical tortuosity, Hydraulic 
permeability is a good indicator of hydraulic tortuosity, and there- 
fore improved relationships between porosity and Formation Factor 


would take into account hydraulic permeability. 
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mt? and Wyble!” and Glen te investigated the effects 
of applied pressure on electrical resistivity values of reservoir 
rocks and sandstone, These rocks of low porosity were placed under 
very high pressures ranging from 1,000 to 12,000 p.s.i.g. The ef- 
fect of pressure on electrical resistivity seems to be much greater 
than the effect on porosity, since electrical resistivity was 
increased much more than porosity decreased. Similar work on sedi- 


ments appears to be nonexistent, 


Mathematical Models 

When selecting mathematical models which best fit a porous 
medium such as a saturated sediment, it is necessary to be aware of 
the influence of the charge on clay particles, clay type, salinity 
of the continuous phase, sediment particle size, shape, distribution, 
and structure, In order to further refine models, awareness of other 
factors such as hydraulic permeability and consolidation pressure 
may be useful, In addition the assumptions of theoretical models 
need special consideration, and extrapolation beyond the range of 
experimentation (in empirical models) is not recommended, 

The electrical conductivity of dispersions has been studied 
extensively. Some studies date back over a hundred years, Only 
those theories that are most likely to apply to sediments are con- 
sidered below, 

їп 1881 ОЛ conducted one of the earliest theoretical 
investigations and derived an equation for the assemblage of spheres 


in a continuous media, The assemblage was a suspension of spheres 
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that were not in physical contact with each other. Maxwell's 


equation, in non-dimensional conductivity terms, follows: 


3 RK¿ 7 2 пка +2 п 
En м Кап - n4 3 ` (16) 








If the spheres are assumed to be non-conductive, i.e, Ky approaches 
zero, Maxwell's equation reduces to: 
2 п 
КИ Ол, (17) 
or in resistivity terms: 
3 - п 
FF = Rn BE (18) 


This equation can be applied to a sediment mixture consisting ofa 

fluid of relatively low electrical resistivity and a network of soil 
grains, not in contact with each other, that are insulators of high 
electrical resistivity. For instance soil granules composed of the 
following common minerals shown in Table 2 would have a high resis- 


tivity. 


TABLE 2 


ELECTRICAL RESISTIVITY OF SOME COMMON SOIL MINERALS 


Mineral Resistivity (ohm meter) 
Clay minerals Tes - ПШ 
Quartz oe - 101% 

7 12 


Саісісе 10 - 10 
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When these soil granules of high resistivity (т) are placed in a 
fluid such as sea water which has a relatively low resistivity (r) 
(0.25 ohm meters), the non-dimensional resistivity of the dispersed 


phase (soil granule) can be expressed as: 


4 
c 
or the inverse relationship: 
S 
Kg - жұт» 75-27% 0. (20) 


ume and Woodside and ios developed the Series and 
Parallel Model Equations which are considered the limiting relation- 
ships for porosity determination by means of conductivity. A 
complete development of these equations is given in Appendix A, 

In the Parallel Equation, the interstitial water and solid 
Soil particle are assumed to be parallel to current flow, The re- 


sulting equation in non-dimensional form follows: 


(Parallel) К. = Ka + n(1 - К), (21) 


The Series Equation evolves when the solid soil particles 


and interstitial water are assumed to be perpendicular to the flow 


of current: 


а 


(Ѕетіеѕ) К = T+ a(K,~ 1)" | (22) 
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Wieder? recognized that other models would fall between the Series 
and Parallel Equations. He proposed a linear combination of the two 
models, Later this linear combination was referred to as the 
Three-Resistor Model by Woodside and ТЕЕ The Three-Resistor 


Model is given below: 


Ko = p (Parallel Model) + (1 - p) (Series Model). (23) 


In this combination p represents the fraction of the sediment 
described by the Parallel Model, 

Woodside and Monero” also proposed a Geometric Equation 
in which the sediment is assumed to be the geometric mean of the 


conductivities of the components of the sediment: 


_ , (len) 
K = К, e (24) 
However if it is assumed that the solid soil particle is non- 


conductive, then K, approaches zero. Consequently the Series 


d 
Equation and the Geometric Equation have no meaning, and the 


Parallel Equation is reduced to: 


Қола, (25) 


and the Three-Resistor Model is reduced to: 


Ra = рп» (26) 


An extensive review of the work up until 1959 including 


investigations by Rayleigh and Runge was conducted by Meredither 
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Не concluded that effects of particle size, shape, distribution and 
structure cannot be ignored if accurate theoretical predictions of 
dispersions are to be made, Meredith developed the Distribution 
Model for treating a concentrated dispersion of spherical and non- 
ЕШ! рагСісіев. The relationships in the Distribution Model 


that are most applicable to a sediment environment follow: 


nonconducting 
‚ _ 8n (n + 1) 
spheres: к= GC и) (УЕДЕ (27) 
nonconducting 
randomly oriented 9a (0521) 
rods: K = (ыз п Оу ; (28) 
rods perpendicularly 
oriented to field 
current: K = n nM : (29) 
m 3-n 


Although the above relationships are presented in non-dimensional 
conductivity terms, it is simple to convert them to Formation Factor 


by using the relationship: 


FF = o (7) 


К | 
3 


Formation Factor-porosity relationships proposed by Maxwell, Wiener, 
Meredith, Woodside and Messmer are presented in Figure 3. Wherever 
Ky is used in formulation it is assumed to be 0,01, Since this value 
is assumed to be very small, 0.01 was arbitrarily chosen. 


Electrical resistivity methods became established as a Geo- 
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Figure 3. 
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Formation Factor-Porosity Relationships Proposed by 
Maxwell, Wiener, Meredith, Woodside and Messmer. 

(In the Series and Geometric cquations the value of 
Қа = 0.01 is assumed; in all other equations Ка = 0.) 
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physical exploration technique in 1927 when Conrad Schumberger 
showed that electrical resistivity measurements of an oil well 
(then called electrical coring), could distinguish between produc- 
tive formation rocks and those that are non-productive, This 
development stimulated extensive research on electrical resistivity 
measurements as an indicator of physical and chemical properties, 
In 1941 КОЛОС defined formation resistivity factor i.e., 
Formation Factor, as: 


FF = - (30) 


ri e 
ЖЕ 


мһеге En is the resistivity of the saturated sediment and г is the 
resistivity of interstitial water. Archie also empirically deter- 


mined an equation relating FF to the formation porosity: 
РЕ = пе, (31) 


where n is the porosity of the formation, and m is the slope of the 
Line in a log-log plot of Formation Factor versus porosity diagram, 
From his experimental work, Archie determined m to be in the range 
of 1.8 to 2.0 for consolidated sandstones, and about 1.3 for loosely 
consolidated sands, 

insane et ШЕР investigated FF to n plots and suggested the 


relation should be: 


FF Yan |! ; (32) 
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where (a) is not always unity, but is often a factor accounting for 
the amount of cementation between mineral grains. ‘Therefore (a) is 
referred to as the cementation factor. From equation (32) the 

"Humble Formula" was developed and is probably the most widely used 


equation in the well logging industry: 


FF 2 0.620 ^' D, (33) 


A comprehensive review of geophysical well logging through 1961 has 
been compiled by Dakhınov. He summarized the factors that affect 
electrical resistivity of a porous media by a general formula, If 


the equation is applied to a saturated sediment, it becomes: 
га = £,(c) £,(n) £,(8) £,(1) (9) Е (к) £,(x,) (35) 


where m is the resistivity of the sediment, 

c the amount of clay/silt in the sediment, 

n the porosity of the sediment, 

S the partial saturation of the sediment, 

Т. the temperature of the sediment, 

Q the cation exchange capacity of the soil mineral, 

ty the resistivity of the soil mineral, 

re the resistivity of the interstitial water, 
Since some of the variables are sometimes omitted, the complexity of 
the above equation is often reduced. For example if the soil miner- 


al (r4) is nonconductive, this term is omitted, In a saturated 


sediment (S) is a constant equal to one, If the environment is 
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saline the interstitial water reduces the effect of cation exchange 
capacity (Q) and cation exchange capacity is usually ignored, How- 
ever if the interstitial water has a low ionic concentration this 
term must be considered, The temperature (T) can usually be measured 
and a correction to electrical resistivity can be made based on the 


following formula; 


есі Пс (1---2577), (35) 
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Ф o Ф * О o 
where rfjg is the resistivity at 25 Centigrade, 


En the resistivity measured at temperature T, 


te the temperature coefficient of resistivity, which usually 
is about 2.5 percent per degree Centigrade, 
In most situations then, the resistivity of a saturated sediment 
(at a set temperature) will be determined by three parameters if the 
environment is saline: ro» €, and n; and by four parameters if the 
ionic solution of the interstitial water is low: ros €» Q and n, 
Dakhnov also obtained a theoretical relation between Formation 


Factor and porosity of uniform nonconducting spheres in a conducting 


fluid: 


1+0.25(1 - п) 3 


FF = e 
L- (1 - 2/2 


(36) 
sinov = suggested a theoretical equation for cemented rocks 
consisting of grains which have the forms of ellipsoids of revolu- 


tion, The formula for the Formation Factor along the axis of 
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revolution is: 


in n 
f + Vo“ ln (1 - у? 
FF EI. . (37) 


In this equation f is defined as: 


18 
f = И, (38) 


(1 - TO - е) 


2 
- П.в, (39) 


е іѕ the eccentricity of an ellipse, a and b being semimajor and 


where 


(D 
l 


semiminor axis of the ellipsoid. The extreme cases of this equa- 
tion are when f = 0 (nonconducting grains of a flat disc shape) 
and when f = 2 (nonconducting grain would be spherical). This 
theoretical work, although developed for the electrical resistivity 
of rocks, combined two additional important factors, that of parti- 
cle distribution which is termed the sorting coefficient (у), and 
a particle shape factor (f). These two factors would also be 
important in a sediment environment. 

Formation Factor, porosity relationships proposed by Archie, 
Dakhnov and Semenov are present in Figure 4. The Humble Formula 


is also presented in Figure 4, 
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Figure 4. 
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Formation Factor-Porosity Relationships Proposed by 
Archic, Dakhnov, Semenov and the Humble Formula. 
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More recent works which are not reviewed previously in the 
literature are discussed below, These recent works are either di- 
rectly applicable or actually conducted on sediments, 

Atkins and mens: have experimentally found m, which they 


call shape factor exponents, for various mineral particles: 


Mineral Particle m 

Na Montmorillonite 3.28 
Ca Montmorillonite 2.70 
1111 6е ІЯ 
Kaolinite 1.87 
Sand 1,60 


They also concluded that (a) the cementation factor in formula 32 
was not unity but was a factor of proportionality which accounts for 
amount of clay to sand ratio, Only if a system consisted of one 
type of mineral particle would (a) be unity. 

It is interesting to note that if m could be found by experi- 
mentation, then for that type of clay only the determination of the 
Formation Factor would be necessary and the porosity could be esti- 
mated, Although this study was conducted in order to understand 
how the clay mineral fraction influences reservoir rocks, the results 
could be very useful in the prediction of porosity from Formation 
Factors in sediments, 

Boyce” took 43 electrical resistivity measurements on sedi- 


ments taken from cores obtained from seven locations in the Bering 


Sea, some of which were 600 nautical miles from each other. Using 
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a consolidometer he squeezed the interstitial water out of these 
sediments and measured the electrical resistivity of this water in 
order to obtain Formation Factors. The porosity of these measure- 
ments ranged oe 58.3 to 87.4 percent. By the method of least 
Squares for an exponential function the following formula was 


derived: 


FF = ОИЕ ой (40) 


Boyce found that this formula had an error of + 15 percent, Не 
stated that this rather large error could be caused by different 
types and percentages of electrical conductors, clay minerals, and 
pore textures. The large error suggests that 600 miles is too large 
an area to attempt to define with one formula, 

Kermabon et TN studied the electrical resistivity of deep 
sediments from 21 long cores taken in a 4 by 10 square mile area in 
the Tyrrhenian Sea, The porosity of these sediments ranged from 
50 to 87 percent, By assuming the resistivity of the interstitial 
water was the same as the bottom sea water, Formation Factors oe 
determined. More than 2500 measurements of porosity and Formation 
Factors were compared, A third degree polynomial curve was found 


to best fit their data relating porosity and Formation Factor 


n(%) = -5.9021 FF? + 40.0416 FF“ - 105.3889 FF + 171.2504 , (41) 


or, 


oa OS 


-= 1.4615 (42) 
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According to Kermabon et al, the above equation fit the data 
obtained with only a + 2.5 percent porosity prediction error, 
Sens studied electrical resistivity measurements of North 
Atlantic deep-sea cores and in situ resistivity measurements on the 
continental shelf around Wales, England. The resistivity of the 
interstitial water was either estimated from bottom water salinity 
or actual core water salinity measurements were taken in the labora- 
EX. A best fit straight line through 99 measurements resulted in 


the following equation: 


ЕЕ = Weenie t, (43) 


Smith noted that the data clearly grouped into two classes: for 


silts and clays of porosities greater than 60 percent, 


FF =n, (44) 
and for sands and coarse silts of porosities less than 60 percent 


palos. (45) 


He concluded that there was no need to use a polynomial forinula of 
the type described by Kermabon et JT but that once a Formation 
Factor is measured an approximate value of porosity could be ob- 
tained from a general empirical equation, and that after selective 
sampling a better value for porosity could be obtained by a more 
refined equation, 


Porosity versus Formation Factor equations are plotted in 
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Figure 5 from Boyce, Kermabon et al., and Smith. Figure 6 shows 
Atkins and Smith's formulations using equation (32) for sodium 
montmorillonite (m = 3.28), illite (m = 2.11), kaolinite (1.87) 
and sand (m = 1,60), 

PESO LAS was able to correlate several sedimentary layers by 
taking electrical resistivity measurements on sediment cores ob- 
tained some distance apart. He did this work both in a lake and in 
a marine environment. Although no formulation between electrical 
resistivity and sediment properties was given, the ability to match 
horizontal layers may indicate a similarity between physical and 
chemical properties with electrical properties of sediments. 

Bouma et En and Chmelik et SEES developed laboratory and 
in situ electrical resistivity equipment for marine sediments. 
Their objective was to make correlations on lithology and geotech- 
nical properties, Although no formulations were given, they at- 
tempted to correlate electrical properties with the chemical, phys- 
ical and engineering characteristics of sediments, Properties such 
as pH, Eh, water content, density, carbonate content, grain size 
analysis, X-ray radiographs, photographs, cone penetrameter and 
vane shear measurements were compared with electrical resistivity 
measurements, Certain trends were apparent from their data.  Elec- 
trical resistivity was indirectly proportional to the water content 
and percentage of clay size present, 

Bani? studied the literature in order to determine the 
feasibility of using electrical conductivity techniques for sediment 


porosity prediction. He concluded that the following points need 





FORMATION FACTOR 


Figure 5. 
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Formation Factor-Porosity Relationships Proposed hy 
Boyco, Kermabon, et. al., and Smith. 
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Figure 6. Formation Factor-Porosity Relationships Proposed hy 
Atkins and Smith for SOdium Montmorillonite, Tllite, 
Kaolinite and Sand. 
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further work both in the laboratory and in the marine environment. 
l, What theoretical or empirical equation best describes 
the electrical conductivity-porosity relationship? 
2. What is the effect of clay on electrical conductivity 
measurements? 
3. How can the conductivity of interstitial water in marine 
sediments be conveniently determined? 
4, What is the best method for restricting the volume over 
which a measurement would be made? 
The first two questions are considered in this section of the paper, 
while the last two questions are investigated in the second section, 
The main motivation behind this literature review and the following 
laboratory experiments is to justify the use of electrical resisti- 
vity techniques as an in situ measuring method to predict porosity 
in a marine environment. Consequently sediment type, interstitial 
water salinity, and formulations applicable to the marine sediment 
environment were studied more extensively than those applicable to 


a fresh water environment, 


EXPERIMENTAL PROCEDURE 


Materials 

The following types of materials were obtained for testing: 

1. Kaolinite clay (Albion Kaolin from Babcock and Wilcox 
Refractories Division), 

2. Illite clay (Grundite from A. P. Green Refractories 
Company), 

3. Montmorillonite clay (Black Hills Bentonite from 
International Minerals and Chemical Corporation), 

4, Providence silt, 

5. Ottawa sand (well-rounded grain shape), 

6. Two types of glacial sands (angular grain shapes), 

7. Four types of marine sediments, 

One type of marine sediment was obtained by a Ewing gravity 
corer and the Deep Ocean Sediment Probe (DOSP) in April 1970 from 
shallow water near Ponce, Puerto Rico. The properties of this sedi- 
ment were extensively investigated by 5. The other three 
types of marine sediments were taken by an Electrical Resistivity 
Measuring System, These samples were taken during the Fall of 1971 
from Narragansett Bay. 

The properties of the marine sediments and the other labora- 
tory materials are presented in Table 3, Co-ordinates of marine 


sediment cores measured in this investigation are given in Table 4, 


397 
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TABLE 3 


MINERALOGY, CLASSIFICATION, GRAIN SIZE, SHAPE, SPECIFIC GRAVITY, 


SEDIMENT CLASSIFICATION 








LIQUID 
NAME OF SEDIMENT MINERALOGY! 7% Sand h% Silt ^ Clay [LIMIT 
E... Kaolinite ши 2 
Quartz 
Illite 












Montmorillonite 
K-Feldspar 
Quartz 






Montmorillonite 


Providence Silt 


Glacial Sand(l- ve 


Glacial Sand (2) 










Quartz 
Plagioclase 





100 









Quartz 










Puerto Rican 
Sediment 


Aragonite 
Calcite 
Quartz 
Plagioclase 


70 - 50 






Narragansett Bay 
(Station A) 


Quartz 
K-Feldspar 
Mica 





41 - 33 


Quartz 
K-Feldspar 
Mica 


Narragansett Bay 


(Station B) 2 





Narragansett Bay 
(Station C) d 


о 

ч 

Co 
1 
t 





l. From X-ray Diffraction analysis (see Appendix С), 
2. 8 = - log, (diameter in millimeters), 
850 = diameter for which 50% of the sample is finer grained. 
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LIQUID LIMIT, PLASTIC LIMIT, pH, AND ORGANIC CONTENT OF SEDIMENTS 























GRAIN SIZE ORGANIC 
: - „ | SPECIFIC PLASTIC CONTENT 
650 6, SHAPE GRAVITY LIMIT pH (%) 
8.38 | 1.81 |platy 2.60 4,1 - 
1 
5-10 
4-5 


3. 6, (Dispersion measure) = 1/2 (¢ 
4. obtained by electron microscope 
and conventional microscope. 


84 ^ 816) 
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TABLE 4 


LATITUDE AND LONGITUDE OF MARINE SEDIMENTS 


I. Puerto Rico Cruise, April 1970 
le 17957.4!' N - 66°37.4' W 
2. 17°57.2' N - 66°27.7' W 
3. 17956.5' N ~ 66%36.6' w 
4. 17056.2' М - 66935.0' и 


5. 17954.4' N - 66935.5' И 


II, Narragansett Bay, Fall 1971 


Station A - 41 34.0' N 71°25.4' W 


71°24.7' н 


Station B =~ 41°29.5' N 


Station C = 41°22.7' N 71930.6' w 
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The details of the tests used to obtain these properties are given 
in Appendix C. These include X-ray analysis, grain size analysis, 
hydrometer, specific gravity, liquid and plastic limit determina- 
tion, pH, and organic content. 

The illite clay, kaolinite clay and Providence silt were 
made as nearly homoionic to sodium ion as possible by washing 
batches of clay in concentrated (1N) solutions of sodium chloride, 
The excess salt was then removed by leaching clays and silt with 
distilled water and the conductivity of the solution was checked 
after each dilution by means of a Beckman Conductivity Bridge 
(Model RC-16B2). When the solution conductivity was that of normal 
tap water, the clays and silt were subsequently dried in an oven at 
105°C, lightly pulverized, and stored in sealed jars. This prepa- 
ration of samples was not conducted on Black-Hills Bentonite since 
it is a sodium montmorillonite clay and therefore is already homo- 
ionic to the sodium ion. The sands were sieved to insure proper 
grain size and then washed in tap water to remove as much ionic 


contamination as possible. The sands were then dried and stored. 


Laboratory Apparatus 
The laboratory testing apparatus was developed in the follow- 
ing stages: 
l. a simple apparatus to determine the feasibility of the 
me thod, 
2. a modified design to obtain variation of key parameters. 
3. redesigned equipment to facilitate testing and allow an 


increased number of tests. 
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Simple apparatus 


The first resistivity cell was constructed from a seven inch 
length of acrylic tubing (I.D. 1438 1ммн ре О О ел5) (GU 
Figure 7). The aluminum top and bottom electrodes were machined down 
to fit the inner diameter of the tubing with O-ring seals to allow 
for a tight seal. A small hole was drilled in the top electrode to 
provide an escape for excess water while the cell was being closed, 
The four side platinum electrodes were installed (see Figure 8) by 
drilling circumaxially two small holes, 1/8 inch apart, and looping 
0.025 in. platinum wire through the holes, thus providing four ap- 
proximate point electrodes spaced 3/4 inch apart vertically along 
the interior of the cell wall, The holes were then made air tight 
by using a GACO sealant, The platinum wire was then soldered to a 
standard copper wire in order to connect the cell to the instruments. 
For testing actual marine sediments the top and bottom electrodes 
were redesigned in order to allow sediment samples of various lengths 
to be used in the cell. (See Figure 9). 

Another simple testing apparatus (see Figure 10) was made by 
placing top and bottom current electrodes on a marine sediment sample 
located in the original core liner used to obtain and store the 
sample. Top and bottom electrodes were machined to the same diameter 
as the core liner, A hole (3/16 in. in diameter) was drilled along 
the perimeter of the top electrode. This hole allowed an insulated 
1/8 inch diameter stainless steel tube to pass vertically through the 
top electrode into the sediment. Located at the end of the stain- 


less steel tubing were two vertically separated, platinum potential 
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The Simple Cell Apparatus. 


Figure 7. 
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Figure 8. The Four Side Platinum Electrodes of the 
Simple Cell Apparatus. 








Figura 9, 


Simple Cell Apparatus Used to Test e 
Of Various Lengths. 
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Testing Apparatus and Electronic Equipment Used 


Figure 10. 


Sediments in the Laboratory. 


to Test Marine 
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electrodes. The wiring from the platinum potential electrodes 
passed through the tubing to the instrumentation, 

Modified apparatus 

A modified design was required in order to test various 
porosities of the fine grained sediments such as silts and clays. 
This design also made it possible to see if consolidation pressure 
altered electrical resistivity measurements to such an extent as to 
mislead the porosity prediction. 

Sedimentation tubes were fabricated with a piston loading 
rod which applied pressure to the sediment, This pressure could be 
applied incremently by adding weights to the loading rod. Each 
tube was constructed from a nine inch length of acrylic tubing 
(I.D. = 1.42 in., O.D., = 1.75 in.). The tubes were bored to insure a 
concentric inner diameter, and allow the piston head and rod assembly 
to slide freely inside the tube. Small holes were drilled into the 
piston head to allow fluid drainage. The piston head was made of the 
same material as the tubes to prevent any binding due to temperature 
changes. A bottom cap of the same design and dimensions as the 
piston head was placed. in the bottom of the tubes to prevent sediment 
loss from the bottom when the pressure was applied. In addition 
micro-pore filters were used above the bottom cap and below the piston 
head to prevent the sediment from oozing out of the fluid drainage 
holes. A PVC (poly-vinyl-chloride) plastic piston rod was screwed 
into the piston head and supported a PVC loading plate. Five side 
Ы кы TN were installed in the same manner described for 


the first cell. Ribbon wire was used to connect the four cells to 
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the instrumentation (see Figures 11 апа 12). By connecting individual 
cells and by using a simple switching device, different electrode 
combinations could be used. The four cells conveniently fit into a 
beaker chamber and were held upright by a plastic support frame. 

Four of these chambers were fabricated at the URI Engineering Instru- 
ment Shop. Details of the design are given in Appendix B. 

The electronic equipment (see Figure 13) used with the simple 
apparatus and modified apparatus was obtained from New London Labora- 
tory Naval Underwater Systems Center, New London, Connecticut. The 
l ma, 1000 Hz signal sent to the current electrodes was supplied by 
a constant current generator and power supply. The voltage differ- 
ence obtained from the potential electrodes was sent to an AC/DC 
converter. The potential difference between electrodes was measured 
from a digital frequency counter and the same signal monitored on an 
oscilloscope. The electronic instrumentation was calibrated and the 
accuracy of the system was determined to be within + 1 percent, 

Redesigned apparatus 

The Е онаи was redesigned in order to reduce the 
time required to consolidate the fine grained sediments. A modified 
version of a pressure vessel designed and fabricated by о: в was 
incorporated. The pressure vessel was constructed from a 3% inch 
length of aluminum pipe which was capped on both ends with cast iron 
pipe caps. The air pressure is regulated and enters the vessel through 
a top fitting. The porous disc sits on the bottom end cap. This end 
cap Nr acts as a reservoir for interstitial water draining out of 


the sediment when the pressure is applied, Pressures up to 19 psig 
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One Cell of the Modified Test Apparatus. 


Figure 11. 
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One Chamber of the Modified Test Apparatus. 


Figure 12, 
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Electronic Equipment Used with the Simple 


and Modified Testing Apparatus. 


Figure 13. 
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can be applied without bubbling air through the porous disc. The 
cell was constructed of thin walled plastic tubing (I.D. = 1.38 in., 
O.D. = 1.50in.) and 5.5 inches long. Four platinum wire 

(dia. = 0.025 in.) electrode rings were cemented in grooved indenta- 
tions on the inside tube wall with epoxy. The vertical separation 
between rings was 0.375 in. The cell was cemented to the porous 
disc with epoxy. The platinum ring electrodes and thermistor (used 
to measure temperature changes) were connected to wiring that fed 
through an air tight seal located at the top of the pressure vessel, 
This wiring connected the cell to the instrumentation, Figure 14 
illustrates the cell, pressure vessel, regulator and gauge assembly. 
Appendix B contains other details. of design, 

New electronic instrumentation allowed continuous electrical 
conductivity measurements during the consolidation process and 
recorded the results. This electronic instrumentation consisted of: 
(1) conductivity signal conditioner, (2) impedance matching amplifier, 
and (3) slow speed single channel recorder. The conductivity signal 
conditioner was manufactured by Honeywell Inc. (Model description 
552022-1002-100-DOO) as part of a water quality measuring system 
(see Figure 10). Since the signal conditioner is basically a solid 
state null balancing A.C. amplifier only the sensor (electrodes) had 
to be modified to allow this system to be used on sediments. The 
signal conditioner supplies a 60 Hz 0.6 ma current to drive the 
sensor artay, and compares the resultant output of the sensor voltage 
electrodes with a stable reference voltage. A differential voltage is 


used to regulate the drive current and it produces an output linearly 








Figure 14. 


Cell, Pressure Vessels, Regulator and Guage 
Assembly Used with the Redesigned Equipment. 
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proportional to the conductivity of the measured medium. This out- 
put voltage can be read on the signal conditioner panel meter, a 
digital voltmeter, or sent to an impedance matching amplifier and 
then to a chart recorder where a continuous record can be maintained. 

Since sediment conductivity varies with operating temperature, 
compensation is required to refer the measurement to standard condi- 
tions (usually 25°с). The thermistor, located inside the pressure 
vessel and attached to the plastic cell, directly alters sensor cur- 
rent, thereby compensating for temperature changes during testing. 

According to Honeywell operating specifications, the calibrated 
accuracy of the signal conditioner and sensor is + 0.50 percent for 
the full scale range. Calibration curves, schematic diagrams and 


detailed operating principles are presented in Appendix B. 


Methods of Testing 
Simple apparatus 


In the simple apparatus first described, an electrical cur- 
rent was allowed to flow between the top and bottom electrodes of the 
resistivity test cell, and a flow pattern developed that was analo- 
gous to Darcy's Law for one dimensional fluid flow, In this study 
when a voltage potential was measured between two platinum potential 
electrodes, the resistivity of the sample was calculated as follows: 


. ABE , , AE 
si laa i el i oe 


where? 
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s the resistivity of the sample, 
A = cross sectional area of the cell, 
a = the voltage electrode spacing, 


AE = voltage difference between voltage electrodes, 
I = the current, 


K 


the geometric cell constant, 
If the current was allowed to flow between the top and bottom 
platinum electrodes, a different flow pattern developed. The sec- 
ond part of Equation 46 is still applicable if the voltage potential 
is measured between the two inner platinum electrodes. The geometric 
cell constant (K) can no longer be determined simply by dividing the 
cross sectional area by the distance between potential electrodes 
but is usually determined empirically by using solutions of known 
resistivities. 

The cell constant was determined by making measurements with 
the cell filled with standard solution of sodium or potassium chloride, 
for which the resistivity may be found in chemical tables, In this 
study, cell н for all the geometrical electrode configurations 
were obtained using 0.01N KCl, 0.1N NaCl, 0.5N NaCl, and 1N NaCl solu- 
tions. All cell constants of simple geometric configuration were 
compared with measured values and less than 1 percent deviations р 
encountered between the two values, 

Since variation of interstitial water was desired in the labora- 
tory both to simulate ocean as well as estuary conditions and to 
observe the effect solution concentration has on the electrical resis- 


tivity of the sediment, sodium chloride solutions of 3,5, 7.0, 17.5 
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and 35.0 parts per thousand were used in the early stages of testing. 
In later testing standard sodium chloride solutions of 0.1, 0.2, 0.3, 
0.5 and 0.6 normality were used as the interstitial fluid of the 
sediment matrix. 

In early laboratory testing using the ӨТІ. се11 m electri- 
cal resistivity of sands was tested using the following procedure: 

The interstitial salt solution was first poured into the cell 
and the sand was added slowly. The sand settled out rapidly in a 
loose packed configuration, The temperature and length of the sand- 
water system was recorded and the top electrode of the cell was then 
secured. The resistance measurements were taken using at least three 
separate configurations of electrodes, Ten measurements of each con- 
figuration were taken and the average value was multiplied by its ap- 
propriate geometric cell constant to give the sediment resistivity. 
Electrode configuration resistivity measurements were compared with 
each other and if variations greater than 3 percent existed, the 
test was re-run. On completion of the first series of tests the 
cell was tapped approximately twenty times and the sand densified 
into a closely packed configuration. Sediment resistivity measure- 
ments were again taken and the length of the sand-water system was 
recorded. On completion of these tests, the top electrode was re- 
moved and the temperature again recorded. Temperature variations 
were seldom more than a degree Centigrade, since the sand and solu- 
tion were usually at room temperature prior to mixing. The sand 
matrix was then dried at 105?C for 24 hours and the weight of sand 


was'obtained for porosity measurements, 
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Clays and silt were tested in a similar manner in the 
simple cell apparatus, with the exception that the clays and silt 
were allowed to equilibrate 24 hours in the respective salt solu- 
tion under a partial vacuum to remove entrapped air. Slurries were 
then added to the cell and allowed to settle under a partial vacuum 
in order to avoid formation of air bubbles. The aa period 
varied depending on the individual sediment, Approximately six or 
seven slurry replenishments with total time periods of 24 to 36 hours 
for settling were required to make a cell ready for testing, Due to 
the long settling periods the clay measurements were only made for 
7.0, 17.5 and 35.0 parts per thousand sodium chloride solution. In 
addition densification of the highly porous clays and silt were not 
performed since tapping the cell caused the resuspension of the sedi- 
ment. The procedure previously explained for sand samples was used 
to measure electrical resistivity and porosity of clays and silt. 

The total time to run sand tests was approximately twenty minutes. 
Although the actual testing time for clays and silt was ten minutes, 
a test cell was occupied for an average period of 24 hours due to 
the long settling periods. 

A natural marine sediment wae studied in the early stages of 
this investigation in order to (a) determine a sediment resistivity 
profile versus porosity and water content, (b) determine the Forma- 
tion Factor by measuring or estimating the resistivity of interstitial 
water, and (c) gain experience for more extensive future testing. As 
previously mentioned, the first marine sediment was obtained by both 


a gravity corer and the DOSP, Certain sections of these cores had 
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previously been segmented into four inch sections for storage, 
Sixteen of these short segmented sections were tested in the 

acrylic resistivity cell, 1 3/8 inch diameter samples were taken, 
The electro-osmosis sample cutter was used in order to reduce sample 
disturbance and core shortening, The samples were transferred to a 
plastic tube of the same dimensions and sealed with rubber stoppers, 
This tube was then placed in a constant temperature bath (25°C) for 
24 hours prior to testing. 

In order to determine water content a sample was taken from 
the excess sediment, Approximately the same amount of sample was 
placed in a test tube and diluted with 100 ml of distilled water, 
The tube was sealed and the mixture was shaken vigorously. After 
the slurry in the tube settled, a conductivity measurement of the 
solution was taken using the Beckman conductivity meter, If it is 
assumed that the number of ions in solution remained constant and 
conductivity changes are linear for small dilution changes, it ís 
possible to approximate the original conductivity of the interstitial 
water, In addition, when a small quantity of interstitial water was 
available, an American Optical Corp. salinity T/C refractometer 
(Model 10402) was used to measure salinity. This permitted the 
determination of the resistivity of the interstitial water, 

After 24 hours the marine sediment sample was removed from 
the temperature bath and placed in the testing cell. Three resisti- 
vity measurements were taken for each sample. These measurements 
corresponded to a segment of the sample which was dissected upon 


removal of the sample from the testing cell. The individual porosity 
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was then obtained for each sample. 

Those marine samples not dissected into small segments were 
analyzed directly in their original plastic corer liner. The stain- 
less steel tube containing the potential electrodes was guided 
through the hole on the top current electrode and vertically pushed 
down through a core sample. Electrical resistivity versus depth for 
the entire length of the core was obtained. Since the tubing passes 
along the plastic corer liner interior wall, negligible sample dis- 
turbance resulted in the central area of the core samples. This al- 
lowed electrical resistivity measurements to be obtained without 
destroying the sediment sample for other engineering soil tests. 

Modified apparatus 

The procedure used to test the effect of consolidation pres- 
sure on clays and silts with the modified apparatus was as follows: 
A е: of Ottawa sand was placed in the four liter beaker chamber, 
This layer allowed the load applied by the loaded cells to be equal- 
ly distributed over the base of the beaker. The support frame con- 
taining the four cells was placed on the layer of sand. A micro- 
pore filter was placed in the bottom of each cell. The chamber was 
then filled with a solution of known normality. A slurry was pre- 
pared in the manner previously described, with a solution of the 
same normality as the chamber solution. Approximately 50 ml of the 
slurry was introduced to the cells. After the slurry settled, a new 
charge was added. This process was continued until all four celis 
in the chamber were filled. A micro-pore filter was placed on top 


of the sediment in each cell and the piston and loading rod was 


62 


placed on the sediment giving an initial pressure of about 
0.01 WM 

The initial vertical length between the loading plate and the 
support frame was measured, The weight of the loadins plate de- 
creased the length, After two or three days the movement stopped. 
The change in the length indicated sediment consolidation, The load 
was then doubled and the procedure repeated, The loads used re- 
sulted in pressures of 0,018, 0.036, 0.073, 0.110, 0.146, 0.292, 
0.366, 0.438 and 0.511 о. Electrical resistivity measurements 
and changes in lengths were taken daily, 

Four to six weeks were required to complete a chamber test. 
During this time period the test apparatus was placed under a clear 
plastic hood in order to prevent evaporation which would change the 
salinity of the chamber water (see Figure 13). A foam rubber mat, 
saturated with water, was placed under the test chamber. A humidity 


dial was placed inside the hood. A thermometer was located in the 


beaker chamber in order to allow temperature corrections to be made, 
At the end of the tests the cells were lifted from the chamber and 
the sediment height in each cell was recorded. All of the sediment 
was then extruded from the cell cylinder and analyzed for solid soil 
weight, water content and porosity. 

Redesigned apparatus 

For the redesigned equipment a calibration curve was obtained 
by filling the cell with standard solutions of sodium chloride of 
known conductivity and measuring the signal conditioner output volt- 


age, Sodium chloride of solutions of 0.1, 0.15, 0.2, 0.5 and 0.6 
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ous were used to obtain calibration curves. Flow patterns, 
cell constants and calibration curves are given in Appendix B. 

The method used to test clays, silt and marine sediments in 
the redesigned test apparatus was as follows: а slurry was prepared 
in the manner previously described. Water from the slurry was ex- 
tracted and placed in the cell to determine an interstitial water 
conductivity, The water was removed from the cell, the slurry was 
poured into the cell and its height was recorded, The slurry was 
analyzed for water content. The vessel was then assembled and pres- 
sure was applied. All sediments were tested at four different pres- 
sures; 0.25, 1.00, 4.00, and 16.00 psi. Each pressurized test ran 
approximately 12 hours or until consolidation at the pressure tested 
had ceased, The pressure vessel was disassembled. The sample height 
and the volume of water in the bottom reservoir was measured, 

The sample was extruded and analyzed for water content and 
porosity. The electrical conductivity of the sediment sample was 
continuously recorded during the entire test period. Sample calcu- 


lations for porosity and water content are presented in Appendix C. 





RESULTS AND DISCUSSION 


The conductivity of kaolinite clay, illite clay, Providence 
silt, Ottawa sand (a) and angular glacial sand (l-a) were plotted 
against the conductivity of the interstitial water (see Figure 15), 
The results from Kermabon et ae: for a clay having a porosity of 
50 percent is also presented in this figure. For each of the sedi- 
ments a linear relation was plotted by the method of least squares. 
The slope of TM lines is FF, the Formation Factor. The porosity 
for each sediment is given on Figure 15. The variation of porosity 
for each type of sediment was within + 2 percent, and this is probably 
the cause of deviation from the least square line in the larger grained 
sediments, In the clays this deviation c also be influenced by 
the cation exchange capacity, The effect that ions on the clay sur- 
face have on the Formation Factor is illustrated by the results 
presented in Figure 16, In this figure the Formation Factor for the 
three clays used in laboratory experiments (kaolinite, illite, and 
montmorillonite) were ЖОО against increasing interstitial water 
conductivity. Though porosity was held constant the curves show 
that the Formation Factors of sodium montmorillonite and illite clays 
increased with increasing interstitial water conductivity from 
1,1 to 3.0 mho/meter. From this point on the Formation Factor be- 
came relatively constant, This result differs from the result found 


by Atkins and Smith; in their work the cation exchange capacity of 
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Figure 15. 


Electrical Conductivity of the Sediment Versus 
the Conductivity of the Interstitial Water for 
Tllite clay, Kaolinite clay, Providence silt, 
Ottawa sand, and Glacial sand. 
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Figure 16. Formation Factor Versus the Conductivity of the 
Interstitial Water for Kaolinite, Illite and 
Montmorillonite Clay. 
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sodium montmorillonite was of no influence on resistivity at an 
interstitial water salinity of 1.1 mho/meter,. Although the 
kaolinite clay data presented in Figure 16 had more scatter, it is 
apparent that the Formation Factor for this clay remained fairly 
constant over the range of interstitial water conductivities used 

in these tests (1.1 to 5.6 mho/meter). These results agree with the 
conclusion made by T that the effect of cation exchange capacity 
on Formation Factor was negligible for kaolinite clay above inter- 
stitial water conductivities of 1.25 mho/meter. 

Sheeler et ino stated that the effect of cation exchange 
capacity on Formation Factor persisted up to interstitial water con- 
ductivities of approximately 6 mho/meter for kaolinite clay, and 
11,8 mho/meter for montmorillonite clay. However examination of the 
data on kaolinite clay presented by Sheeler et al, indicates that 
above an interstitial water conductivity of 1.6 mho/meter the rate of 
increase of the Formation Factor is so small that it could be consid- 
ered constant. In the case of the montmorillonite clay the low 
porosities could account for the increase in Formation Factors with 
interstitial water conductivities up to 11.8 mho/meter. The poro- 
sities of clays used by Sheeler et al. ranged from 48 to 67 percent, 
while the porosities used in this investigation ranged from 
87 to 94 percent, 

Although Boyce” offered cation exchange capacity as a pos- 
sible explanation of the variation in his data (in which the lowest 
interstitial water conductivity encountered in the cores was ap- 


proximately 5 mho/meter), the results of the present investigation 
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indicate that above 3 mho/meter the effect of cation exchange capa- 
city on Formation Factor is negligible for sediments, Five mho/meter 
is equivalent to a salinity of 32.6 parts per thousand at 25 degrees 
Centigrade, while 3 mho/meter is equivalent to 18.5 parts per thou- 
sand, Although no cation exchange capacity values were given by 
Kermabon et ал and mu after conducting laboratory experiments 
on marine clays having porosities of approximately 50 percent, both 
concluded that there was no appreciable change in Formation Factor 
when using interstitial water conductivities ranging from 1,4 to 

13 mho/meter, 

Another interesting result of the present investigation is 
that below 1,5 mho/meter the Formation Factor for montmorillonite 
clay was less than one, This means the sediment has become more 
conductive than the interstitial water. This was also noted by 
Sheeler et шне? when working with montmorillonite clay of lower po- 
rosities, This fact would be extremely important when working 
with montmorillonite clay with low interstitial water conductivities, 
since a correction factor would have to be developed before Formation 
Factor could be considered an accurate index of porosity, 

The Formation Factor and the porosity data of a specific 
sediment are given in Appendix C, A relationship using Equation 32 
was developed, If Equation 32 is re-written the following relation- 


ship is obtained; 


Log FF = log a - m log n. (32) 
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From the above relationship it is evident that (m) is actually the 
slope of the line and (a) is the value of the Formation Factor at a 
porosity of 100 percent or the intersection of the y-axis of a 
log-log plot. Using a linear regression technique these coefficients 
(m and a) were determined by the least squares method. An example 
of a complete data scatter plot relating Formation Factor to poros- 
ity for kaolinite clay is presented in Figure 17, The dashed lines 
in Figure 17 are the + 2 percent error lines located from the line 

of best fit. 

Data for all the sediments tested in this investigation are 
plotted as Formation Factor versus porosity in Figure 18. The lines 
of best fit are drawn for each sediment tested; the sands, silts, 
clays and marine sediments. The length of the line depends on the 
maximum and minimum values tested, This figure represents a graphic- 
al method of predicting porosity from Formation Factors obtained by 
electrical resistivity measurements, The ranges of porosity and 


Formation Factor obtained are presented in Table 5, 


TABLE 5 


THE POROSITY RANGE AND FORMATION FACTOR RANGE OF THE 
SANDS, SILT, CLAYS, AND MARINE SEDIMENTS TESTED 


SEDIMENT TYPE POROSITY RANGE (%) FORMATION FACTOR 


Sands 
Silt 
Clays 


Marine sediments 
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Figure 18. Formation Factor-Porosity Relationship 
for all Sediments Tested. 
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The following data for each sediment tested is presented in 
Table 6: number of data points used in the statistical analysis, 
slope, intercept, standard error of estimate, percentage of data 
points that fall within + 2 percent error lines, and the maximum 
and minimum Formation Factor and porosity values, The techniques 
used in this statistical analysis are discussed in Appendix C. 

Because of the ranges of Formation Factor and porosity (see 
Figure 18) the discussion will be facilitated by dividing the 
comparisons into the following three groups: 

a) sands and silt, 

b) clays and silt, 


c) clays, silt and marine sediments, 


Sands and Silt 

Formation Factor versus porosity is plotted for the five sands 
and Providence silt in Figure 19, Formation Factor for the sands and 
silt in this figure ranged from 3,12 to 5.90, while porosity ranged 
- from 26 to 48 —3 The Formation Factors for Ottawa sands and 
glacial sands (l-a) and (l-b) were very similar at the different poros- 
ities, while the more naturally occurring sediments (Providence silt, 
glacial sand (2) and Narragansett Bay sand (Station C)) were each 
uniquely different. Narragansett Bay sand, although actually a marine 
sediment, was considered under sands and silts since 97 percent of the 
particles were in the sand size range. 

The lines of best fit for glacial sand (2) and Providence silt 


at porosities of 41 to 46 percent visually suggest that electrical 
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resistivity decreases with increasing particle size, Providence silt 
which has a mean diameter of 0.016 mm had a higher Formation Factor 
than glacial sand (2) which has mean diameter of 0.50 mm. At a 
porosity of 43 percent, the Formation Factors were 4.9 апа 3,7 
respectively. The distributions of these two sediments were similar, 
As discussed later, the shape of the silt was similar to the sand, 
that is angular. in addition to particle size, structure was also 
variable in these two sediments. At 43 percent porosity Providence 
silt is densely packed and glacial sand (2) is loosely packed. How- 
ever the results discussed later show that structure of sands does 
not influence electrical resistivity measurements independent of 
porosity. Consequently it appears that particle size did account 

for the difference in the Formation Factor of glacial sand (2) and 
Providence silt. 

These results are in agreement with Meredith: who reported 
that mixtures of particles that were smaller in size by 10 diameters 
resulted in slightly higher Formation Factors than mixtures of larger 
particles, The diameters of glass beads used in Meredith's work 
ranged from 6 mm to 0.02 mm while the size of sand was not given, 

The data of Wyllie and pcs agreed that the effect of different 
sized glass spheres on Formation Factor was noticeable only when the 
difference in particle size was large. The diameters of the glass 
spheres used in their study ranged from 3 mm to 0,033 mm. 

An additional size test was conducted in order to see if small 
differences in particle size of a sediment could be detected.  Elec- 


trical resistivity measurements of Ottawa sand (a) with a diameter of 
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0,70 mm and Ottawa sand (b) with a diameter of 0.42 mm resulted in 
Formation Factors of 3.6 and 3.5 respectively. The porosities, 
particle shapes, and standard deviations of these two sands were 
approximately the same. The work of Meredith" on glass beads agrees 
with this work on sediments that the effect of particle size on 
conductivity is negligible when the size difference is small, 

The effect of distribution on electrical resistivity can be 
seen in Figure 19 by comparing Narragansett Bay sand with glacial 
sand (l-a). The Formation Factor of Narragansett Bay sand which had 
a larger distribution was higher than the Formation Factor of glacial 
sand (l-a). For example at a porosity of 40 percent the Formation 
Factors were 4,36 and 3,85 respectively. The standard deviation of 
the particle size for Narragansett Bay sand was 1.10 while glacial 
sand (l-a) was 0.34, Particle size and shape were approximately the 
same (see Table 3). These results agree with the work of a 
on sedimentary rocks, and the work of Meredith? on glass beads, 
Meredith found that Formation Factor was higher for a mixture of 
spheres 6 mm and 0.75 mm in diameter than for spheres 6 mm in dia- 
meter. 

At lower porosities the effect of distribution was found to 
be even more pronounced. At a porosity of 32 percent the Formation 
Factor of Narragansett Bay sand was 5.70 ий sand (l-a) was 
4.77. Therefore the difference between the two was 0.51 at 
40 percent porosity and increased to 0.93 at 32 percent porosity. 

The increased assortment of particle sizes causes the current to take 


a more tortuous path and thus the sediment is more resistive, Аз 
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the sediment becomes densified (lower porosity), the variety of 
particle sizes has an even greater effect on the tortuosity. 

The Narragansett Bay sand had a small quantity of platy mica 
particles which could be thought to account for some of the dif- 
ference between Narragansett Bay sand and glacial sand (l-a). How- 
ever as discussed later in this paper, particle shape differences in 
sediments do not seem to significantly influence electrical resisti- 
vity measurements, 

Figure 20 a to f are photographs of Ottawa sand, glacial 
sands l-a, l-b, 2, Narragansett Bay sand and Providence silt respec- 
tively. An evaluation of angular and spherical shapes as an influ- 
ence on electrical resistivity can be made by comparing Ottawa sand 
with glacial sand (l-a). Although Ottawa sand is spherical in shape 
and the glacial sand (l-a) is angular in shape, the Formation Factor 
versus porosity relationship is very similar, At a porosity of 
35 percent Ottawa sand had a Formation Factor of 4,20 while glacial 
sand (l-a) had 4,35. Size, distribution and structure were approxi- 
mately the same for these two sands. On Figure 19 95 percent of the 
data points lie within + 2 percent of the line of best fit, Conse- 
quently the difference between the line of best fit for Ottawa sand 
and glacial sand (l-a) is not significant. These results indicate 
that electrical resistivity is quite similar for angular and spheri- 
cal shapes. 

Wyllie and Gregory" obtained a Formation Factor of 4.2 for 
glass spheres and spherical Ottawa sand, 4.4 for angular beach and 


creek sands, and 4.7 for lucite triangular prisms at 35 percent 





Figure 20-a. 





E O, mm. — 


Ottawa Sand. 
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Figure 20-b. 





Glacial Sand (1-a). 
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Figure 20-c. 





Glacial Sand (lb). 
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Figure 20-d. Glacial Sand (2). 
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Figure 20-e. 





Narragansett Bay Sand. 
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Figure 20-£. 





Providence Silt. 
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porosity. The sphericity was 1,0, 0.8, and 0,6 respectively. The 
Formation Factor in their work seems to increase as sphericity de- 
creases. However the particle size was not held constant, The 
particle size of the glass spheres was 3 mm in diameter and ee 
angular prisms were 0.25 mm in diameter, As discussed previously 
this difference in particle size is large enough to influence the 
resistivity. The Formation Factor would be higher for the triangu- 
lar prisms since the diameter was 12 times smaller than the glass 
spheres, 

The influence of structure on the electrical resistivity of 
sands was tested by making resistivity measurements before and after 
densification. The electrical resistivity and porosity were changed 
by densification in the same proportion, The Formation Factors of 
the five sands at 40 percent and 32 percent porosities, the percent- 
age change in the Formation Factor and the measure of spread are 
Piven in Table 7, For a 20 percent decrease in porosity the Forma- 
tion Factor increased 19 to 24 percent. As shown in Table 7, the 
Formation Factors of the two sand types with a large distribution 
changed a little more than the sands with a small spread in particle 
size. 

Tortuosity factors as defined by Equation 13 are given in 
Table 8 for loosely and densely packed sands. Tortuosity increases 
as density increases for those sands that have a large variety of 
particle sizes. Densification does not seem to affect the tortuosity 


of those sands with a small spread in particle size, 





TABLE 7 


THE FORMATION FACTORS AT 40 PERCENT AND 32 PERCENT POROSITY, 
THE PERCENTAGE CHANGE IN FORMATION FACTOR AND THE 
MEASURE OF SPREAD FOR FIVE SANDS. 








(1) (3) 

M FF FF % FF (2) 6, 
40% 32% inercased 

Ottawa sand 3.73 4.60 19 0.29 
Glacial sand (1-а) 3„85 4.77 19 0.34 
Glacial sand (1-Ъ) 3.73 24.60 19 0.50 
Glacial sand (2) 3.96 5.00 21 1.31 
Narragansett Bay 4,35 5.70 24 1.10” 


sand (Station С) 





(1) FF 0% - Formation Factor at 40% porosity. 


4 
FF. - FF 
(2) % FF increased = __32% 40% (100). 
FF 424, 


(3) Sy (Dispersion measure) = 2 (9:4 - AME where §,, and $e are 
cumulative percentages of the sample which is coarser grained than 84 
and 16 percentiles respectively. The f scale is defined as -log, 


(diameter in millimeters). 


# 
* 
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TABLE 8 
COEFFICIENTS OP PERMEABILITY, FORMATION FACTORS, POROSITY, TORTUOSITY, 
AND THE PRODUCT OF THE FORMATION FACTOR AND THE COEFFICIENT OF PERMEA= 
BILITY FOR FIVE LOOSELY AND DENSELY PACKED SANDS. 


Name of Sediment Coofficient Formation Porosity Tortuosity K FE 


of Permea- Factor in) (wi — VER 
bility K, EF). 
(em/sec) 
Ottawa Sand 0.239 3.95 42.0 12202 0.848 
(loose) 
Ottawa Sand 0.154 5.30 2745 1.20 0.816 
(dense) 
Glacial Sand 0,0891 3.43 45.0 1.24 0.306 
(1-a) (loose) 
Glacial Sand 0.0502 5.00 30.2 17023 0.251 
(1-a) (dense) 
Glacial Sand 0,0224 3.32 44.0 1921 0.0744 
(1-b) (loose) 
Glacial Sand 0.00699 4.80 31.2 1.22 0.0336 
(1-b) (dense) 
Glacial Sand 0.115 3.50 44.0 1.24 0.403 
(2) (loose) 
Glacial Sand 0,0413 5.85 27.2 1.26 0.252 
(2) (dense) 
Narragansett Bay 0.0879 4.30 40.5 1.32 0.378 
Sand (loose) 
Narragansett Bay 0.0248 5.75 31.8 1.35 0.143 


Sand (dense) 
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Coefficients of permeability (Kp) and the corresponding 
Formation Factor and porosity are also presented in Table 8. The 
coefficient of permeability as obtained by a variable head perme- 
ability test, was reduced as the sands densified and Formation 
Factors increased. The coefficient of permeability decreased as 
the mean particle size decreased. It decreased more for densified 
sands with a large distribution of particle sizes. As the porosity 
of a sand decreased, the coefficient of permeability decreased 
logarithmically and the Formation Factor increased logarithmically. 
Due to the cancelling effect of these two logarithmic variations, 
the product of Formation Factor times coefficient of permeability of 
densely and loosely packed sand can be plotted on a linear scale. 
This linear assumption, along with a Formation Factor versus poros- 
ity plot allow the coefficient of permeability to be obtained by 
measuring only the Formation Factor, Thus by knowing any one of the 
three variables (FF, n, Kp), the other two variables can be predicted 
from standard curves, The product of the Formation Factor and the 
coefficient of permeability (FFKp) is presented for all the sands on 
Table 8. The relationship between porosity and FFKp for loosely and 
densely packed sands is plotted in Figure 21. By assuming a poros- 
ity of 36 percent and using Figures 19 and 21, the FF FFKp and the 
Kp for the five sands used in these tests are given in Table 9, 

The sands that are more common in nature have the greater slopes 
in Figure 21. For these sands a small change in porosity results in a 
larger change in FFKp than for the sorted sands, Consequently these 


graphs predicting porosity and permeability from electrical resisti- 


89 


POROSITY (%) 


46 


44 


42 


40 


38 


34 


32 


30 


29 


OTTAWA SAND 

GLACIAL SAND (lo) 
GLACIAL SAND (Ib) 
GLACIAL SAND (2) 


х + Ро 


l 2 P 4 5 6 ft ‚В 
(FOROMATION FACTOR) (PERMEABILITY) FFK (¢™/sec) 


NARRAGANSETT BAY (STATION C) 


Porosity Versus the Product of the Formation Factor and 


Figure 21. 


the Coefficient of Permeability (FFKp) for Five Sands. 


| 





90 


TABLE 9 


FORMATION FACTORS, COEFFICIENTS OF PERMEABILITY AND THE PRODUCT OF 
FORMATION FACTOR AND THE COEFFICIENT OF PERMEABILITY AT 36 PERCENT 
POROSITY FOR FIVE SANDS. 





Name FP 363 FFKp 36% Кр 36% 
(cm/sec) (cm/sec) 
Ottawa Sand 4.12 .834 „202 
Glacial sand (l-a) 4.26 .274 ‚ 0643 
Glacial Sand (1-b) 4.12 . 048 .0117 
Glacial Sand (2) 4.40 .321 .0730 


Narragansett Bay Sand 4.95 .258 „0521 
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vity measurements would be very useful in a natural environment, 


Clays and Silt 


A Formation Factor-porosity diagram with curves of best fit 
for kaolinite, illite, and sodium montmorillonite clays and 
Providence silt is shown in Figure 22, This figure includes over 
200 measurements made with all three versions of the laboratory test 
equipment, The statistical data for this figure is presented in 
Table 6. 

The results of this investigation on clays and silts suggest 
that particle shape does not influence electrical resistivity inde- 
pendent of porosity. Kaolinite clay (0.6N NaCl) and Providence silt 
have almost identical Formation Factors over a porosity range of 
51 to 61 percent although the former is platy shaped and the latter 
is angular in shape. These shapes, enlarged 2000 and 10,000 times 
by an electron microscope, can be seen in Figures 23-a, 23-b, 20-f 
and 24, The particle distribution for the two sediments was similar; 
the measure of spread for kaolinite was 1.81 while Providence silt 
was 1.21. The particle size of Providence silt was approximately 
five times larger than kaolinite. It was suggested in the discussion 
on results of sands that a particle size difference of about 10 times 
is necessary before size difference is dere, electrical re- 
sistivity measurements, In the results and discussion on cation ex- 
change capacity it is indicated that cation exchange capacity would 
not be a factor influencing the electrical resistivity measurements 


of these two sediments since the salinity was high (0.6 normal sodium 
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Figure 22. Formation Factor-Porosity Relationships for 
Kaolinite, Illite and Montmorillonite Clays 
and Providence Silt. 
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Figure 23. 
a,b 





Photographs of (a) Kaolinite Clay Enlarged 2000 
Times by Electron Microscope, (b) Kaolinite Clay 
Enlarged 10,000 Times. 


93 





ee 


Figure 24, 





Photograph of Providence Silt Enlarged 10,000 
Times by Electron Microscope. 
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chloride). Since both would have a flocculant structure at this 
salinity, structure would not be a variable. Since shape is the 
only major variable between the two, and since the Formation Factors 
are nearly identical in this porosity range, it appears that parti- 
cle shapes such as platy or angular do not influence electrical re- 
sistivity. The results for illite clay further substantiate this 
conclusion, Although illite clay is of a different mineral content 
from kaolinite, the Formation Factors were again almost identical. 
This result indicates that this difference in mineral content has 
little if any effect on electrical resistivity measurements,  Illite 
clay has a platy shape that can be seen in Figure 25-a and 25-b. 
The mean particle size of illite clay is about the same as for 
kaolinite clay (five times smaller than Providence silt). The Forma- 
tion Factor of illite clay and Providence silt is nearly identical 
though shape differs and appears to be the only major variable. The 
opposite is true of kaolinite clay and montmorillonite clay. Although 
these two have a very similar shape, the Formation Factor versus 
porosity relationship for these two is very different (Figure 26-a, 
and 26-b). For example at a porosity of 83 percent, Formation Fac- 
tor of kaolinite is 1.45 while montmorillonite is 2.00. As discussed 
later this difference is probably due to particle size and distribu- 
tion. The results for illite and montmorillonite clays reinforce the 
conclusion that the influence of particle shape on electrical resis- 
tivity measurements of sediments appears to be small, 

The standard deviation of particle size for illite clay was 


2.99 and might appear to be a variable when compared with Providence 


! 








Figure 25. Photographs of (a) Illite Clay Enlarged 2000 Times 


a,b by Electron Microscope, (b) Illite Clay Enlarged 
10,000 Times. | 
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Figure 26. 
a,b 





Photographs of (a) Montmorillonite Clay Enlarged 
2000 Times by Electron Microscope, (b) Montmoril- 
lonite Clay Enlarged 10,000 Times. 
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silt. Although this standard deviation is twice that of Providence 
silt, this difference is probably not significant. As a descriptive 
measure of spread, standard deviation is less accurate when the 
sediment grains are very Eine, This is due to the assumption, in 
the hydrometer method of analysis, that all particles are spherical 
and do not aggregate while settling. This assumption can be used 
with reasonable accuracy for most TRITT 0.005 mm in size, 
However below 0.005 mm particles vary widely in shape and aggregate 
more readily. The formula for the standard deviation of the parti- 
cle size 64 is based on the 84th percentile diameter (see Table 3), 
Eighty four percent of the Providence silt particles were coarser 
than 0.006 mm while 84 percent of the illite clay particles were 
coarser than 0.0006 mm. It is apparent that error could have been 
introduced by aggregation and non-spherical particles in the illite 
clay analysis, 

Descriptive measures such as phi median diamter (850) or 
phi deviation measure (6 4) determined by a standard hydrometer meth- 
od of analysis cannot be used for montmorillonite clay. Approximate- 
ly 80 percent of the particle sizes are below 0.005 mm in this clay. 
Other methods of analysis such as optical equipment or centrifugal 
methods are required. However time and cost make these methods im- 
practical. Spread or Зее нонок can be estimated from the per- 
centages of sand, silt and clay that make up the sediment. Mont- 
morillonite clay which has a large percentage of clay sized parti- 
cles could be expected to have a large spread in particle sizes. A 


large distribution would be expected to increase Formation Factor. 
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Montmorillonite was more resistive than illite and kaolinite, al- 
though size probably influenced this result also. As discussed 
earlier large‘ distributions increased Formation Factor for sands 
when mean particle size was constant. Just as was found in the 
results on sands and silts, particle size appears to influence elec- 
trical resistivity measurement of clays although porosity is constant, 
As stated earlier higher Formation Factors were obtained for mont- 
morillonite than for kaolinite. One of the two variables (in addi- 
tion to.distribution) is particle size. The mean particle diameter 
of montmorillonite is reported to be<< 0,001 mm while kaolinite is 
0.003 mm. According to os the particle diameters for mont- 
morillonite clays ranged from 0.0016 to 0.00005 mm, These measure- 
ments were made using a super-centrifuge operating at various speeds 
up to 50,000 rpm. This very large size difference in particle size 
most likely accounts for some of the difference in the Formation 
Factor of montmorillonite and kaolinite, At 82 percent porosity, 

the Formation Factor of these two clays was 2.10 and 1.48 respec- 
tively. It is doubtful that the mineral difference of these two 
clays influenced the Formation Factors, since both are non-conductive 
and the interstitial water was 0.6 normal sodium chloride. 

Smaller differences in mean particle size (i.e. Providence 
silt (0.016 mm) and kaolinite (0.003 mm)) did not appear to influ- 
ence the electrical resistivity measurements. As seen in n 22 
the Formation Factors were nearly identical over the porosity range 
of 51-61 percent for these sediments. Because of the larger sized 


particles, Providence silt was not in a stable condition above 





100 


61 percent porosity. If the particle size is smaller (as in 
kaolinite clay) the sediment exists in a stable condition at poros- 
ities up to 84 percent. Electrical resistivity measurements were 
conducted up to a 91 percent porosity by producing a slurry. The 
data of Atkins and smith’? for kaolinite and sodium montmorillonite 
clay is also presented in Figure 22. The Formation Factor-porosity 
data agrees well with the results obtained in this investigation. 
However in their work the difference between kaolinite and montmorillo- 
nite was attributed to shape and values for a given shape factor (m) 
were determined. As previously discussed the results of the present 
study on clay sediments suggest that the difference between the two 
is due to particle size and distribution. 

Meredith! after studying the electrical conductivity of 
dispersions (spheres and discs), concluded that errors as high as 
10 to 20 percent were due to shape differences. Semenov found that 
disc shaped particles in rock resulted in higher Formation Factors 
than rocks of spherical particles at equal porosities. The results 
of the presemt. study indicate that shape is not a factor influencing 
the electrical resistivity of sediments, The earlier discussion on 
sands and silts showed that differences between spherical and angu- 
lar particles were not detected by electrical resistivity measure- 
ments, The results on clays and silts indicato that nearly identi- 
cal electrical resistivity measurements resulted from sediments with 
angular or platy particle shapes as the only apparent variable, 

The effect of structure on electrical resistivity can be seen 


in Figure 22 if kaolinite (0.6N NaCl) is compared with kaolinite 
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using 0.2 molar sodium tetra phosphate solution (Quadrafos) as the 
interstitial water of the sediment, The structure naturally en- 
countered in a saline environment, a flocculant structure (sce 
Figure 2), resulted in a lower Formation Factor than when Quadrafos, 
a common deflocculating agent, was used to obtain a dispersed struc- 
ture. The difference between the two increased as porosity de- 
creased, For example at a porosity of 80 percent the Formation Fac- 
tor for kaolinite (0.6N NaCl) was 1.33 while kaolinite (Quadrafos) 
was 1.55. At a porosity of 52 percent the Formation Factors were 
2.30 and 3.60 respectively. 

In the dispersed structure the clay particles are in a 
preferred orientation, that is one particle tends to be oriented to 
the next particle in the same manner, However in the flocculant 
structure the clay particles are in a randomly oriented pattern to- 
ward each other. According to Brindley?" preferred orientation will 
bring out or increase the intensity of X-ray diffraction patterns of 
M spacing (d-spacing) reflections. The d-spacing is the spacing 
between crystal lattice planes. 

A side test was conducted to check if the particle orienta- 
tion assumed previously could be validated by X-ray diffraction 
techniques. Two X-ray diffraction slides of equal thickness were 
prepared under identical conditions for kaolinite clay with 
0.6N sodium chloride as the interstitial water and kaolinite clay 
with 0.02M sodium tetra phosphate. 

Using the General Electric diffractometer (model XRD-5S) in 


the Metallurgical Research Laboratory at the University of Rhode 
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Island, each sample was scanned from 8 to 26 degrees. The proce- 
dure used is presented in Appendix C. This scan covered the two 
major d-spacing values of kaolinite clay; 7.2 and 3.6 Angstroms, 

The X-ray diffractograms of the two samples are presented in 

Figure 27, The vertical axis is a logarithmic indication of inten- 
sity corresponding to an approximate quantitative index of the 
mineral present, Although many factors affect intensity of the dif- 
fracted X-ray beams, orientation of particle and the reinforcement 
of intensity caused by preferred orientation or the reduction of in- 
tensity caused by random orientation of the same mineral particles 
would be a major factor determining the magnitude of intensity. The 
difference in intensity between the two samples is illustrated in - 
Figure 27, An average intensity count was made at the two major 


d-spacing values, The results are shown in Table 10. 


TABLE 10 
THE AVERAGE INTENSITY COUNTS OF THE TWO MAJOR D-SPACING 
VALUES OF TWO STRUCTURES OF KAOLINITE CLAY 


Sample name 7.2 d-spacing 3.6 d-spacing 


Intensity (counts) | Intensity (counts) 





Kaolinite (0,6N NaCl) 8758 8884 


Kaolinite (Quadrafos) 47333 47842 


The results show that the kaolinite clay (Quadrafos) reinforced the 
intensity when compared with kaolinite clay (0.6N NaCl). These re- 


sults indicate that kaolinite clay (Quadrafos) was of a dispersed 
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structure while kaolinite clay (0.6N NaCl) was of a flocculant 
structure, 

The salinity of marine and estuarine environments results in 
clay deposits formed by flocculation, eee high lateral or 
vertical pressure or perhaps even disturbance could cause the floc- 
culant structure to collapse to a dispersed structure. These results 
Suggest that the electrical resistivity measurements would change as 
the structure changed from flocculated to dispersed. 

The effect of consolidation pressure on the conductivity of 
clay and silt sediments can be seen in Figure 28, These curves il- 
lustrate that the electrical conductivity is inversely proportional 
to the applied pressure in the range from 0.25 to 16.00 psi. In this 
range the Formation Factor increased in proportion to porosity de- | 
crease, Consequently consolidation pressure seems to have no effect 
on electrical resistivity of sediments independent of porosity. 

Figure 28 presents four curves showing illite clay consoli- 
dating at pressures of 0.25, 1.00, 4.00, and 16.00 psi and the re- 
sulting conductivity values of 2.7, 2.5, 2.0, and 1.4 mho/meter 
respectively. Similar results obtained for Providence silt, kaolinite 
clay and montmorillonite clay, can be found in Appendix С, Porosity 
was determined at the beginning and end of the consolidation ргеѕ- 
sure tests, 


As discussed in the Review of the Literature several studies 
18,19,20 


D 


on formation rocks have reported that Formation Factor in- 
creases much faster than the decrease in porosity. The pressures 


used (1,000 to 12,000 psig) were much greater than the 16 psi used 
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in а study. It would seem possible that sediments would be more 
easily influenced by overburden pressure than rocks. However the 
greater increase of Formation Factor over porosity decrease was not 
P nna on sediments at the low consolidation pressures used in this 
Brody. 

A typical conductivity versus time recording obtained by 
the redesigned equipment is presented in Figure 29 for kaolinite 
clay (0.6N NaCl). Most tests were run for 12 hours. At the times 
designated as ti and t, in the figure, a Formation Factor was ob- 
tained and compared with porosity values determined on the initial 
slurry and after consolidation. These two values were compared with 
Formation Factor-porosity values for kaolinite (0.6N NaCl) from 
Figure 22, At ty and t, the Formation Factors were 1,64 and 3,27 
respectively. These values fall within the + 2 percent porosity 
error lines for kaolinite (0.6N NaCl) in Figure 22, The values for 
conductivity of the sediment sample were changed to Formation Factor 
by dividing the conductivity of the interstitial water by the con- 
ductivity of the sediment. By comparing these Formation Factors 
with those of Figure 22, values of porosity were obtained and used 
on the vertical scale of Figure 29. Since porosity and void 


ratio (e) are uniquely related by: 


the void ratio was included as an addition on the vertical axis. 


For saturated sediments a relationship exists between the void ratio 
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Figure 29. Conductivity, Formation Factor, Porosity, Void 
Ratio and Water Content Versus Time of Kaolinite 
Clay at a Consolidation Pressure of 16 psig. 
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and water content (w) that can be obtained by knowíng the specific 


gravity of the sediment: 
wea. (48) 


The specific gravity (G) of kaolinite clay is 2.60. Water content 
was calculated and also included on the vertical axis of Figure 29. 
Therefore Figure 29 is a continuous record of the change in Forma- 
tion Factor, porosity, void ratio and water content of a kaolinite 
clay sediment as it is consolidated by pressure of 16.00 psig over a 
twelve hour period, The total test time, consolidation pressure, 

and for times ti and с, the Formation Factors, porosities, void 
ratios, and water contents for all sediments tested by the redesigned 
equipment are presented in Appendix C. 

The continuous conductivity readings by the redesigned equip- 
ment can be extremely useful, If accurate Formation Factor-porosity 
relationships are determined for a sediment as they were in this 
study, then from the continuous conductivity readings (Formation Fac- 
tor) the porosity, void ratio, and water content can be obtained at 
any point during the monitoring. Using the void ratios obtained at 
the end of a test, a void ratio versus the log of the consolidation 
pressure is plotted in Figure 30. In classical soil mechanics this 
type of curve is very useful in predicting consolidation, If a 
Formation Factor-porosity curve of a sediment is sufficiently accu- 
rate, then the redesigned equipment could be used as a consolido- 


meter, Curves such as Figure 30 could be obtained by simply contin- 
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Formation Factor, Porosity and Void Ratio Versus 
the Log of the Consolidation Pressure of Illite 
Clay. 
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ually recording the electrical conductivity of the sample and 
incrementally increasing the pressure after set periods of time 
without removing the sample. Then by using the Formation Factor- 
porosity relationship the porosity and consequently void ratio 
could be obtained and plotted versus the corresponding pressure, 
Currently tests are being conducted with a more advanced teflon 
cell in the redesigned equipment, This cell reduces side wall 
friction and allows more extensive and elaborate testing to be 
conducted. 

Since the water driven out of the sample during consolida- 
tion was collected in the reservoir of the redesigned equipment, it 
was possible to obtain an approximate value of the coefficient of 


permeability (K ). Thus by using Darcy's Laws 


t 
м 5 (49) 


where 

K_ is the coefficient of permeability (cm/sec), 

q the quantity of water collected in the reservoir during 
the test em, 


t the total time of the test (sec), 


e? 


the average length of the sample (obtained by measuring 
the initial length of sample and adding it to the final 
sample length, then dividing by two) (cm), = 

Ah the pressure used in the test (cm), 


A the cross-sectional area of the sample (en 
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An example of the results of this type of calculation is presented 

in Figure 31. In this figure the void ratio at the end of a test is 

plotted against the coefficient of permeability for the illite clay. 

From this curve it is possible to obtain other coefficients of perme- 
ability from Formation Factors, if the Formation Factor-porosity 


relationship has been developed, 


Marine Sediments 

For the four different marine sediments tested 
90 to 100 percent of the Formation Factor-porosity data points fell 
within + 2 percent of the line of best fit. This indicates that 
electrical resistivity measurements can accurately predict porosity 
of marine sediments, The values that fell outside the + 2 ы 
error lines were unusual samples usually containing pieces of coral 
or shells, 

Two vertical profiles of a marine sediment core sample ob- 
tained by using the simple laboratory test apparatus are presented in 
Figures 32-a and 32-b. Porosity, water content, and sediment resis- 
tivity (corrected to a temperature of 25°C) are shown. The points 
on the curves represent discrete measurements, The horizontal dashed 
line indicates where the core was cut to facilitate handling and 
storage, 

The electrical resistivity of this sediment was inversely pro- 
portional to the water content and porosity. An exception is illus- 
trated in Figure 32-а, А large shell approximately two inches in 


diameter was located at a depth of 20 inches, At this location the 
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Figure 31. Formation Factor, Porosity and Void Ratio 
Versus the Coefficient of Permeability of 
Illite Clay ө 
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Two Vertical Depth Profiles of Porosity, Water 
Content and Electrical Resistivity of a Marine 
Sediment. 
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resistivity of the sediment increased sharply, but the water content 
T porosity remained fairly constant, These sharp changes in re- 
sistivity are indicators of structure and fabric changes along the 
core sample.. One 2À inch diameter core showed sharp increases in 
resistivity. Wet sieving a five inch section of this core through a 
40 mesh sieve (the section with the highest electrical resistivity). 
produced over 320 grams of shells, coral and other calcium carbonate 
material. 

Formation Factor-porosity data obtained by using both the 
simple and redesigned equipment for the four marine sediments and 
the lines of best fit are presented i Figure 33, Since some of 
this data was obtained using slurries, the vertical dashed lines in 
Figure 33 indicate the porosity below which the sediments exist in 
a stable condition, 

Discussion of the factors which influence the Formation 
Factor-porosity relationship for патос ОЕ is difficult be- 
cause with the exception of Narragansett Bay sand (Station C), the 
marine sediments tested were mixtures. The mixtures of particles 
found in a marine sediment can be seen in Figure 34. This photo- 
graph (after a 575X electron microscope enlargement) of the sediment 
obtained at Narragansett Bay (Station B) shows a large angular sand 
grain and a variety of plate, rod and angular shaped particles. The 
purpose of this picture is not to determine particle shape which has 
been found to have little effect on electrical resistivity measure- 
ments of sediments, but to illustrate the wide variety of particle 


sizes and shapes found in a marine sediment, Measures such as 850 
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Figure 34. Photograph of a Marine Sediment from Station B 
Enlarged 575 Times by Electron Microscope. 
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and 97 can be obtained for marine sediments but the values for 
different sections of a marine core sample vary so widely that 
little can be concluded from these measurements. Consequently 
sediment classification, that is percentage of sand, silt, and clay, 
is more useful than particle size and shape when observing the ef- 
fect marine sediment particles have on electrical resistivity meas- 
urements, For example it was observed that the Formation Factor of 
sediments with a higher sand content change less with a change in 
porosity than of sediments with a lower sand content, The reduc- 
tion in negative slope of the Formation Factor-porosity curve for 
Narragansett Bay sand (Station C) can be seen in Figure 18, If the 
line of best fit is extrapolated to the vertical axis at a porosity 
of 100 percent, those sediments with a higher sand content resulted 
in higher intercept values. Table 11 illustrates these points by 
showing the Formation Factor intercept values and the negative 
slopes along with the percentages of sand in the marine sediments 


tested, 


TABLE 11 


FORMATION FACTOR INTERCEPT VALUES, NEGATIVE SLOPES, 
AND PERCENTAGE SAND OF FOUR MARINE SEDIMENTS 


(Intercepts) 
Sediment FF at 100% porosity Slopes j|Sand % range 


Narragansett Bay 


(Station C) Lei 
Puerto Rícan 

sediment ұсы 
Narragansett Bay 1 15 

(Station B) pe 
Narragansett Bay 





(Station A) 
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The marine sediments which contained a higher percentage of 
clay were stable at a higher porosity than the sediments wíth less 
clay. This point is illustrated in Table 12, which shows the per- 
centage of clay and the highest porosity at which the marine sedi- 


ments existed in a stable condition, 


TABLE 12 


PERCENTAGE CLAY AND THE HIGHEST POROSITY AT WHICH 
THE FOUR MARINE SEDIMENTS WERE STABLE 


Clay (%) Highest 


Porosity (%) 


Puerto Rican marine sediment 


Narragansett Bay (Station A) 


Narragansett Bay (Station B) 


Narragansett Bay (Station C) 





The sediments that can exist at a higher porosity result in a larger 
range of Formation Factors since the range of porosity at which elec- 
trical resistivity measurements can be made is greater, 

Data such as coefficient of permeability, porosity, void ratio 
and water content versus time for various consolidation pressures is 


presented in Appendix C. 


Model Equations 


The coefficients of the empirical model equation (Equation 32) 
were determined for each sediment. A similar method is used by well 


logging engineers to determine empirical model equations for different 
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rock formations, Table 13 presents the empirical models for each 
sediment and the range of Formation Factors and porosity determined 
in this investigation, 

Since 94 percent of the data points used to obtain the line 
of best fit fell within + 2 percent porosity error lines (the other 
six percent fell within + 4 percent), а оше can be consid- 
ered good empirical prediction equations for each individual sedi- 
ment,  Kermabon et a obtained a porosity error line of 
+ 2,5 percent with 2500 data points used to determine a prediction 
equation from 21 cores taken from a 40 square mile area in the 
Tyrrhenian Sea, The prediction equations developed by 
Kermabon et арам воусе, ^ апа stent? (presented in the Review of 
the Literature section as Equations 41, 40, and 43 respectively), 
were tested with this data and predicted the porosity of Providence 
silt, kaolinite clay (0.6М NaCl), illite clay and Narragansett Bay 
sediment (Station A and B) within + 5 percent in the 
50 to 70 percent porosity range. Although the "Humble Formula", 
Equation 33, is the most widely used prediction equation for rock 
formations, it is not at all applicable to sediments, 

The two factors which influenced Formation Factor-porosity 
equations the most in the saline environment of these laboratory 
experiments, were the mean size of particles (if the mean size dif- 
fered by at least ten diameters) and the distribution or spread of 
particle sizes, The data in Figure 18 is grouped into two distinct 
areas: in one area were the sands, and in the other area were the 


silts and clays. As discussed previously, classification into 
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TABLE 13 


THE COEFFICIENTS OF THE MODEL EQUATION, POROSITY RANGE 
AND FORMATION FACTOR RANGE OF THIRTEEN SEDIMENTS. 


Coefficients of the 





Model Equation? Porosity Formation 
“Name of Sediment FF = an” Range (%) Factor Range 
—0„97 
Ottawa sand ГЕ = 1.55 n 26 — 43 3.5 = 9.7 
-0.96 
Glacial sand (1-а) ГЕ = 1,57 п 30 - 47 3.2 - 5.0 
-1.10 
Glacial sand (1-b) ГЕ = 1.35 n | 31 142 3-427 
-1.10 
Glacial sand (2) FF = 1,45 n 31 - 46 odd - 5.4 
-1.20 
Narragansett Bay (Station C) РР = 1.47 п 32 - 40 4.3 - 5.7 
Kaolinite clay (0.6N NaCl) FF i sod 91 - 92 1.2 - 3.8 
-1.8 
Illite clay FF = 1,08 n : 99 - 81 1.6 - 3.3 
з | | 
Montmorillonite clay FF =n 11 81 - 96 1.4 - 2.1 
-1.24 
Kaolinite clay (Quadrafos) FF = 1,05 n 4} - 94. 1,1 - 3,1 
-1.75 
Providence silt FF = 1,11 n 4] - 62 A49 31542 
-2.22 
Narragansett Bay (Station A) FF =n 47 - 82 1.9 - 5.0 
-1.75 
Narragansett Bay (Station B) ГЕ = 1.15 n 44 - 74 2.0 - 4.9 
-1.77 
Puerto Rico marine sediment FF = 1,36 n 46 — 74 2.1 ~ 5.0 


*On double logarithmic paper (m) of the model equation is the 
slope of the line of best fit and (a) is the intercept of the equation’s 


line with the vertical axis at-:100 percent porosity. 





121 


percentages of sand, silt and clay generally indicate mean particle 


size and spread. A prediction equation such as 
FF = 1,6/n, (50) 


would approximate the porosity of the sands of this study to within 


+ 5 percent porosity. A prediction equation such as 
2 
FF » l.1/n , (51) 


would best approximate the porosity of silts and clays to within 

+ 5 percent porosity with the exception of montmorillonite. Mont- 
morillonite, which has a high clay content (78 percent) resulted in 
increased slope value (m). This probably can be explained by the 
fact that this sediment has a much smaller mean particle size and 
greater distribution than found in sediments containing less than 


70 percent clay. A prediction equation such as: 


РЕ 


Vn, (52) 


or 


FF 


Ma, (53) 


would approximate montmorillonite and similar clays better. Of 
the four marine sediments tested in this study the porosity of 
three was predicted within + 5 percent by the equation which best 
approximated clays and silts, while the porosity of Narragansett 
Bay sand (Station C) was predicted within + 5 percent by the equa- 


tion which best approximated sands. It appears that porosity can 
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be predicted within + 5 percent by obtaining an electrical resisti- 

vity measurement and knowing the major type of sediment present: 

(a) sand (b) silt and clay (c) predominantly clay. Table 13 sug- 

gests that a more accurate prediction is possible by determining 

the electrical resistivity and the porosity of enough samples to 

obtain a Formation Factor-porosity curve (at least ten samples with 

a wide range of porosity), and developing an equation that approxi- 

mates that specific sediment. This is easily accomplished using the 

redesigned equipment. This equation could then be used to predict 

the porosity (within + 2 percent) of the area by measuring only the 

electrical resistivity of the sediment and its interstitial water, 
As discussed in the Review of the Literature section, most 

of the theoretical model equations were developed for dispersions 

or porous mediums such as rocks. However the following observations 

compare the results of this study with the theoretical work in the 


literature, 


1. All data obtained in this investigation lie within the 
Parallel and Series Equation (Equations 21 and 22). 

2. The equations that considered particle shape such as Maxwell, 
Meredith, and Semenov (Equations 18, 27 and 37 respectively) 
agree with the Formation Factor-porosity relationships of this 
study; however the results on particle shape in this investi- 
gation show that in sediments shape is not ӨЙ ИИО factor. 

3. The non-dimensional conductivity term (Ка) is unnecessary in a 


saline environment (above 18.5 parts per thousand at 25°C) since 





4. 


5. 
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this term approaches zero. Equations that consequently have 
no meaning in a saline environment are: Equations 22 and 24. 
However ina fresh water environment the term Кі A 0 апа these 
equations may be applicable, especially if a clay with a high 
cation exchange capacity such as sodium montmorillonite is 
present, 

Kaolinite clay with 0.02M sodium tetra phosphate solution as 
interstitial water followed the Maxwell Equation (Equation 18) 
for non-conducting spheres though the particle shape of this 
clay was platy (see Figure 23-a and 23-b). The structure of 
this sediment was dispersed and therefore the grains were not 
in direct contact with each other. As stated earlier shape does 
not seem to influence the resistivity of sediments. These 
results suggest that the Maxwell Equation is a good prediction 
equation for sediments of dispersed structure regardless of. 
their shape, The dispersed structure of the kaolinite clay is 
close to one of the boundary equations (Parallel Equation). 
This is as expected since the dispersed orientation is an 
extreme case. 

In the general theoretical formula (Equation 34) for resistivi- 
ty of a porous medium such as a sediment, the term f(c) is a 
factor that depends on the amount of clay/silt in the sediment. 
The results of thís investigation indicate that the term f(c) 
is more dependent on the mean particle size and distribution 


or spread than on the clay/silt shape or mineral content, 


La 


2. 


CONCLUSIONS 


Formation Factor-porosity relationships were determined for 
three clays (kaolinite, illite and montmorillonite), Providence 
silt, four sands and four marine sediments, Formation Factors 
ranged from 1.1 to 5,9 while porosities ranged from 

26 to 93 percent. The results show that the porosity of these 
sediments correlated well with electrical resistivity measure- 
ments, In fact 94 percent of the data points fell within 

+ 2 percent of the line of best fit, while the other six percent 
fell within + 4 percent of the line of best fit, 

Mathematical model equations and a graphical presentation of 
the data showed that the data grouped into two classes, 

a) sands and b) silts and clays. In the mathematical equation 
(FF = an ^), the intercept value (a) was higher for sands (1.6) 
than for clays (1). The slope value (m) was approximately 
minus one for sands, minus two for silts and clays, and minus 
three or four for sediments having at least 70 percent clay 
sized particles. This description of (a) and (m) excludes 
kaolinite clay partícles in a dispersed structure, While the 
intercept value was similar to the other clays (1.05), the 
slope (m) was closer to that of sands than clays (-1.24). 
Kaolinite clay particles in a dispersed structure followed the 


theoretical equation of et for non-conducting spheres, 
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The cation exchange capacity influenced the electrical conducti- 
vity measurements of sodium montmorillonite and illite clays 
when the conductivity of the interstitial water was l.l to 

3.0 mho/meter. Above 3.0 mho/meter this influence was negli- 
gible. The cation exchange capacity of kaolinite was negligi- 
ble over the range of interstitial water salinities tested 
(1.1 to 5.6 mho/meter). 

The mineral difference between kaolinite and illite clays did 
not influence the electrical resistivity measurements in this 
study. 

As the particle size of sediments decreased, electrical resis- 


tivity increased. When the mean particle size of a sediment 


was 30 times larger than that of a similar sediment at the same 


porosity, the Formation Factors were considerably different, 
However small differences in particle size (about 5 times) 

were not detectable by electrical resistivity measurements, 

The distribution of particle size was found to affect electri- 
cal resistivity measurements; the Formation Factor was found to 
increase as the spread in particle size increased. 

The particle shape of the sand, silt and clay sediments tested 
did not appear to influence the electrical resistivity measure- 
ments in this study. Spherical shapes were compared with angu- 
lar shaped particles and angular shapes were compared with 
platy shaped particles, 

The structure of sands did not influence the electrical resisti- 


vity measurements independent of porosity. Densely packed sand 
| 


| 
| 
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was compared with loosely packed sand, 

Consolidation pressure (for normally consolidated sediments) 
up to 16 psig did not influence the electrical resistivity 
measurements of the clays and silts tested independent of 
porosity. 

The structure of two kaolinite clay samples did affect electri- 
cal resistivity measurements, especially at lower porosities, 
Гес есед structure had a much higher Formation Factor 
than the dispersed structure at 52 percent porosity. 

The marine sediments varied so much in particle size, shape 

and distribution within core sections that percentage sand, 
silt and clay is a more useful classification. For the marine 
sediments with higher sand content the change in electrical re- 
sistivity was less with a change in porosity than the sediments 
with less sand, 

Unusual marine sediment samples which fell outside the 

+ 2 percent error lines of the Formation Factor-porosity line 
of best fit were found to contain large amounts of shell or 
coral. 

The electrical resistivity measurements made with the re- 
designed laboratory equipment in this study allowed perme- 
ability, tortuosity, and void ratio-log consolidation pressure 
to be determined. Though further study is required to validate 
this as a method, electrical resistivity measurements by the re- 
designed equipment could provide a fast, easy and inexpensive 


method of testing in the field of submarine soil mechanics, 
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Тһе redesigned laboratory equipment which consisted of a cell, 
an electrode array and a signal conditioner is safe, simple, 
easy to operate and economical ($500.00). Also required for 
these tests are a standard volt meter and recording equipment 


(a simple strip chart recorder 15 sufficient). 


II 


AN ELECTRICAL RESISTIVITY MEASURING SYSTEM 
FOR THE DETERMINATION OF POROSITY OF 


MARINE SEDIMENTS IN SITU 





INTRODUCTION 


An Electrical Resistivity Measuring System was designed, 
fabricated, and tested as a method of predicting porosity of in situ 
marine sediments, The system consists of electrode arrays, an elec- 
tronic measuring circuit, and a FM multiplex telemetry data link, 
These three components will be added to the Deep Ocean Sediment 
Probe (DOSP). | 

The DOSP (see Figure 35) was designed and built by the 
University of Rhode Island under contract to the Naval Underwater 


Systems Center (NUSC), Newport-New London пе 


). The original 
mission of the DOSP was to obtain the velocity and attenuation of 
compressional sound waves through the ocean sediment, This was 
accomplished by inserting four probe legs, located at the corners of 
the base structure, incrementally into the sediment to a distance of 
five feet, In one of the legs an acoustic sparker source provides a 
signal which is received by hydrophones located in the other legs. 
In addition, a thin walled corer, centrally located equidistant from 
the probe legs, is used to collect a sediment sample for laboratory 
correlation with in situ acoustic and temperature data obtained at 
various probe penetrations, With a current design depth of 

5,000 feet, the DOSP provides an excellent stable ocean floor plat- | 
form for in situ sediment measurement sensors, In addition to the 
Electrical Resistivity Measuring System, other systems currently be- 


ing added to the DOSP are a Nuclear Gage (Rose and в) for 
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Figure 35. Photograph of the Deep Ocean Sediment Probe 
(DOSP). 
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sediment density measurements, a geochemical probe, and a multi- 
coring unít to obtain longer cores, 

The ability of the DOSP to obtain samples with minimum dis- 
turbance for correlation with acoustical data is a major step for- 
ward, However, certain variables such as; (1) a measure of dis- 
turbance while penetrating, (2) temperature changes and possible 
release of gases in solution while returning the sample to the sur- 
face, and (3) sample handling on deck, enroute to the laboratory, 
and in the laboratory, are not monitored and therefore still leave 
an area for improvement. The Electrical Resistivity Measuring Sys- 
tem will allow continuous resistivity readings to be taken while the 
corer barrel is penetrating the ocean es (both inside and out- 


side the corer barrel) and also while the sample is being returned to 


the surface, It will be possible to estimate the degree of disturb- 
ance during these operations and in addition to predict sediment 
porosity and water content, Resistivity readings could also be taken 
after sectioning (just before storage) and later in the laboratory 
prior to testing analysis, and thus allow a measure of disturbance 


during these handling procedures, 





REVIEW OF LITERATURE 


Most electrical resistivity field equipment was developed by 
geophysicists for oil or other mineral A DA re- 
viewed almost all the geophysical well logging equipment used be- 
fore 1962. A detailed description of all the electronics and 
electrode arrays used is included. When and where each electrode 
array is best used while surveying is also discussed. Books and 
pamphlets of more recent equipment used in geophysical well logging 
are available from Schlumberger Ltd., P. O. Box 2175, Houston, Texas, 
77001. Electrical resistivity measuring equipment and interpreta- 
tion of electrical resistivity logs are included, Pier suggests 
that the unique hardware which has been developed for the radical 
changes in pressure and temperature that occur in oil well bore 
holes could have useful application in the ocean as well, 

The application of two different geophysical methods in order 
to determine the physical properties of land soils has been described 
by Borowczyk and mM S These authors combined electrical 
resistivity and radioactive isotopic methods in one probe called the 
"Universal Probe", This probe consisted of a steel tube with elec- 
trical resistivity electrodes located in the tube head. The probe 
was introduced into the soil by either a vibrating hammer or drill 
rig - pile driver which used a tripod for vertical support. While 
sinking the probe into the ground at intervals of 10 or 20 cm., elec- 
trical resistivity measurements were taken. After these measurements 
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radioactive isotopic transmitting and detecting equipment were low- 
ered down the steel tube (inner diameter 3.5 cm.) and the intensity 
of the scattered gamma radiation was determined. The authors con- 
cluded that the Universal Probe allowed accurate soil water contents 
and porosities to be predicted by nuclear measurements, soil bed and 
thickness identification by electrical resistivity measurements, and 
approximate soil strength by the penetration resistance measurements 
made during probe insertion, These conclusions were arrived at by 
correlating the measurements made by the probe with known strati- 
graphy and soíl properties obtained by borings. 

The electrical resistivity surveying techniques that were 
first developed by me were used in a land survey by Pertét, °° 
With a simple apparatus consisting of batteries, a double commutator, 
a null indicator potentiometer, and the necessary cables, reels and 
steel electrodes, he was able to determine the depth of a marine clay 
which underlies the lower Mississippi Valley. At four check points 
the electrical resistivity data agreed with boring data within 
3.5 percent. By using electrical resistivity surveying techniques, 
‚it was possible to use boring intervals of 5,000 feet instead of the 
usual 1,000 feet. This resulted in an 80 percent reduction in time 
and a 59 percent reduction in cost when compared with surveying the 
site only by boring techniques, 

The standard geophysical laboratory or field electrical re- 
sistivity equipment used on land can be used to measure the resisti- 
vity of marine sediments. The first resistivity measurements made 


on marine sediments were conducted on board ship after cores had been 
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taken. Cell modification was required to handle sediments instead 
of formation rocks. More extensive equipment modifications are 
necessary in order to measure resistivity of marine sediments 

in situ. Since bore holes are not provided, a penetration type 
probe is required, Deep penetration is not usually possible since 
power sufficient to drive a probe deep into the bottom is not avail- 
able on the deep ocean floor. Since a stable platform is usually 
unavailable, distance of probe penetration is difficult to measure 
accurately. Either the electronics and power must be encased in a 
water tight pressure vessel or long lengths of electrical cable must 
be provided. 

Boyce" measured electrical resistivity of Bering Sea cores on 
board ship and in the laboratory. The original plastic core liner 
served as the веш? A copper ring was securely glued on a circular 
jacket. Four silver alloy needles were soldered to the copper ring. 
The needle electrodes were geometrically spaced 90 degrees apart 
around the ring. When the circular jacket was clamped to the core 
liner, the needle electrodes penetrated the core liner and made con- 
tact with the sediment. Four of the circular jackets were vertically 
spaced along the core liner, This configuration is similar to four 
ring electrodes in which the current electrodes are the outer two 
rings and the inner rings are the potential electrodes, The elec- 
tronics consisted of a FRL Resistivity meter (Type 3320, manufactured 
by Scientific Services Laboratories, Inc., Dallas, Texas). This 
instrument was basically a potentiometric nulling circuit, The de- 


vice was nulled manually and the resistance of the sediment was di- 
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rectly read on the meter. The resistance reading was then correc- 
ted for temperature and multiplied by a cell constant to give a re- 
sistivity value corrected to 25° Centigrade, Interstitial water was 
Squeezed out in a press and salinity was measured by an induction 
salinometer, Although sediment resistivities were obtained both on 
board shíp and in the laboratory, the resistivity of the interstitial 
water obtained in the laboratory was used to obtain Formation Factor. 
for both the sediment at sea and in the laboratory. 

The apparatus used by Pautot^? to measure electrical resisti- 
vity of sediment cores taken from both the sea and a lake environ- 
ment was a quadripole electrode arrangement which was moved along 
the vertical axis of the core. The quadripole electrodes were of 
the Wenner arrangement, A Tellohm resistance meter measured the re- 
sistance directly. The resistance reading was multiplied by a geo- 
metric constant to obtain the resistivity. 

aa {п order to determine bottom types (sand, silt 
and clay) designed, built and tested a device that measured the 
in situ Formation Factor at the water sediment interface. The For- 
mation Factor measuring electrodes consisted of a mutually shared 
copper ring electrode placed on the perimeter of a non-conducting 
plastic disc. Two copper point electrodes were placed on each side 
in the center of the disc. The point electrodes were insulated 
from each other. The disc was placed at the water-sediment inter- 
face to locally isolate the water and sediment for separate resis- 
tance measurements. Since the upper and lower electrode geometry 


was identical, the ratio of the resistance at the terminals was 
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equal to the Formation Factor. The electronics of the system 
consisted of a ratio bridge circuit that was manually nulled to 
directly read the Formation Factor. These Formation Factor measure- 
ments are based on the assumption that the resistivity of the water 
immediately above the sediment is the same as the resistivity of 

the interstitial water in the adjacent surface sediments. The 
temperature of the sediment and water are also assumed to be the 
same. 

Since electrical resistivity measurements on cores are time 
consuming and expensive, Bouma et sae and Chmelik et ap s modi- 
fied an electrical well logging system to be used with a specially 
designed in situ probe. The in situ probe consisted of a heavy 
stainless steel tube 12 feet in length with stabilizing fins welded 
to the upper end. The electrodes consisted of brass rings mounted 
to the nose of the stainless steel tube. Rings of PVC insulated the 
electrodes from each other and from the stainless steel rod. The 
electrode wiring which was encased in a water tight cabling was 
brought up through the stainless steel tube to the well logging 
electronics located on board the research vessel, The electrodes 
were designed to operate in the normal mode and in a guard mode. 
The normal mode is an array in which the distance between the measur- 
ing electrode is large compared to the distance between the current 
source and the closest measuring electrode. The guard mode is an 
array in which a thin center current electrode is placed between two 
guard current electrodes, This configuration causes the current flow 


from the center electrode to be limited to a horizontal plane rather 
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than spreading out. The measuring electrodes are usually located 
between the guard electrodes and the center current electrode, 

The in situ probe was lowered on a steel cable and allowed to 
fast fall the last ten feet by releasing the brake, A spoked wheel 
(18 inches in diameter) which contained magnets in the wheel hub was 
leosely fitted to the stainless steel probe tube. As the probe 
penetrated the bottom the wheel remained at the water sediment inter- 
face, Magnetic proximity switches were placed at one foot intervals 
on the probe. Depth of probe penetration was obtained by an over- 
riding spike signal caused by the closing of the proximity switches 
as they passed through the magnetic field located in the wheel hub, 
As the probe was extracted an excitation current and the return sig- 
nal were sent over a separate cable to obtain the resistivity 
measurements, 

Electrical resistivity measurements were made on the cores 
using the same electrical well logging system, but the electrodes 
were mounted to a sliding carriage which was in synchronous movement 
with a recorder. Cores were vertically cut in half and one section 
was placed in a rack. The electrodes, consisting of brass pins, were 
lowered into the sediment, The equipment was energized and the elec- 
trodes raked the sediment as the carriage moved down the length of 
the core section, A continuous electrical resistivity reading versus 
core length was recorded, 

It is difficult to interpret the results obtained with this 
equipment since temperature is not accounted for and interstitial 


water salinity and bottom water salinity were not given, Consequent- 
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ly ras Factors could not be obtained and the work could not 
be compared with other research work, 

New equipment proposed by Bouma et ӘКЕ eliminates the need 
of the wheel depth indexing device, This is accomplished by having 
a series of point electrodes arranged in a spiral on a fiberglass 
barrel so that there is a 1.1 inch vertical spacing between elec- 
trodes. After the probe penetrates the bottom only sequential 
switching would be required to obtain an electrical resistivity versus 
depth measurement, 

Kermabon et іы developed both laboratory and in situ equip- 
ment. The resistivity cell used in the laboratory experiments had an 
electrode array consisting of four electrode rings. A 16 Hz square 
wave alternating current was fed through the two outer ring electrodes 
and the potential difference was read across the two inner ring elec- 
trodes. The same 16 Hz square wave current was used in the in situ 
resistivity probe. The electrode array was a four point electrode 
system mounted on a plastic pad located at the end of the probe. The 
outer two point electrodes were the current electrodes while the in- 
ner two electrodes measure the potential difference, 

The electrical resistivity probe weighed 100 kilograms in 
water and was 13 meters in length. Since the probe was designed for 
deep sea work, all the electronics — encased in a large pressure 
housing. The electronics consisted of a battery power supply, an 
electrical resistivity circuit, a depth monitoring circuit and a mag- 
netic tape recording unit. The operating principles of the in situ 


probe were as follows: 
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l. The probe was lowered to 10 meters from the sea floor. 
2. The probe was released and allowed to free fall and make bottom 
penetration, 
3. Upon penetration the reference platform remained on the surface 
of the sea floor while the probe continued its downward thrust. 
4, After the completed penetration the probe was extracted; elec- 
.trical resistivity measurements were made with the four point 
electrode array and the penetration was measured via a mechan- 
ical pulley hookup to the reference platform, 
5. The electrical resistivity versus penetration data was stored 
on tape for later playback on a X-Y recorder in the laboratory. 
To obtain Formation Factors, Kermabon et sd assumed the electri- 
cal resistivity of interstitial water was the same as bottom water 
resistivity. Bottom water resistivity values were obtained by 
measuring bottom water salinity and temperature, It also had to be 
assumed that the temperature of the bottom water was the same as 
that of the sediment to a depth of 13 meters. This equipment has 
been transferred from SACLANT ASW Research Centre Laboratory at 
La Specia, Italy, to the Marine Science Laboratories University 
College of North Wales for further testing. 
smith’? reported on resistivity measurements made in the 
laboratory and in situ by using commercially available electrical 
prospecting equipment and the Schlumberger 4-electrode arrangement. 
To obtain Formation Factors from the in situ measurements, the re- 
sistivity of the bottom water was measured before penetrating the 


sediment, then this resistivity value was assumed to be the same as 
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the interstitial water. Іп the laboratory the resistivity of the 
interstitial water was calculated from knowledge of its salinity. 

If porosity of marine sediments is to be accurately predicted 
by an Electrical Resistivity Measuring System, a convenient way of 
measuring the resistivity of the interstitial water must be found, 
One method involves the use of a consolidometer to squeeze out the . 
interstitial water in order to determine salinity which is then 
easily converted to resistivity. The squeezing process is laborious, ' 
time-consuming, and requires large sediment samples. Since samples 
must be analyzed in the laboratory at a later date, desiccation may 
occur and introduce error. Also it is possible that squeezing causes 
error by filtering the water, particularly if clay is present, 

A simpler method is to measure the bottom water salinity and 
assume that this is the same as the salinity of the interstitial 
water. The accuracy of each of these methods has been considered, 
Boyce" squeezed the interstitial water from cores and found that 
this.salinity was not the same as the salinity of the bottom water, 
He found that salinity of the interstitial water varied with loca- 
tion and depth. However since the measurements were made three to 
ae months after the cores were taken, Boyce stated some desicca- 
tion probably occurred, The maximum difference between the inter- 
stitial water salinity from that of overlying bottom water was 
2.3 parts per thousand while the typical difference was 1.2 parts 
per thousand, 

Siever et > squeezed the water from twenty-two cores . 


taken in six different areas and concluded that the variations within 
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a core can be significant. However his data indicates that in most 
cases the average value of interstitial water salinity for a core 
was very close to the bottom water salinity for the area. 

Smirnov 1 studied a large amount of empirical data which was 
collected in the 1950's and 1960's (including the data of 
Siever et ЫК О, and concluded that the salinity of the interstitial 
water of marine sediments does correspond to the salinity of the 
bottom sea water. He had previously proved this point theoretical- 
ly. Data by аа agrees that the difference in salinity be- 
tween the two is small, After studying Pacific red clay, Atlantic 
red clay, calcareous ooze, claying mud, varved clay, and lacustrin 
mud, the largest difference between the salinity of the interstitial 
water and that of bottom water was 0.8 parts per thousand. The data 
of Kermabon et йа also suggests the error introduced by assuming 
these salinities to be the same is very small; for 2500 resistivity 
measurements with depths down to 13 meters the Formation Factor- 
porosity line of best fit had only a 2.5 percent error. 

In order to predict the porosity of marine sediments accurate- 
ly by an Electrical Resistivity Measuring System, it is important to 
use an electrode array that restricts the volume of the sediment 
sampled. Past work with in situ measurements has been done with 
electrode arrays that encompass a large volume of sediment, 

Bouma et Ee found that less detail was recorded with the in situ 
probe than on deck measurements with laboratory equipment, The 
M. detail was attributed to the in situ probe's electrode array 


which sampled a larger area than desired. Consequently equipment 





142 


was planned that was expected to reduce the area sampled. 

Kermabon et a? used the four-electrode system by ee 
in the nose of the resistivity probe. A detailed derivation of the 
Wenner electrode configuration is presented in Appendix A, A cur- 
rent was passed between the outer two electrodes, and the differ- 
ence in potential between the inner two electrodes was measured. 
Theoretically, if the medium surrounding the electrodes was of 
uniform, isotropic, homogeneous resistivity, the calculation of the 
true resistivity could be obtained from the knowledge of the current, 
the potential difference, the geometric arrangement, and the separa- 
tion of the electrodes. Therefore if the electrodes were equidistant 
and the distance between electrodes was large in comparison with elec- 


trode penetration into sediment, the true resistivity would be: 





ср = EN, (54) 
where Er is the true resistivity assuming a homogeneous, isotropic, 
uniform media; 
a, the electrode spacing; 
AV, the potential difference between the voltage electrodes; 
I, the current; 
K, the geometric coefficient for the particular electrode 
array. 
However, the medium around the electrode array is not usually 
homogeneous, isotropic, and uniform; therefore the value computed in 


equation 54 becomes the apparent resistivity (г,). The apparent re- 
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sistivity may be larger than, less than, or equal to the true resis- 
tivity. The difference between the apparent resistivity and the 
true resistivity will depend greatly on the volume of sediment 
encompassed by a specific electrode array while the resistivity 
measurement was being obtained, 

Using the skin depth conduction theory, ene showed the 
volume sampled by a Wenner electrode configuration was directly pro- 
portional to the distance between current electrodes and inversely 
И portional to the frequency of the current, and the conductivity 
of the medium, Then, using a frequency of 1000 Hz and a current 
electrode spacing of 3 cm,, Ball determined the volume of material 
through which a 1000 Hz A.C. current passes as a function of the 
conductivity of the material, Higher frequencies are undesirable 
due to A.C. impedance effects. Using this electrode conr КОЕ 
and a normal range of marine sediment conductivities 
(0.015 to 0,035 mho/cm.) the volume of material that the 1000 Hz 
A.C. current could encompass would range from 0.6 to 1.2 cubic ЕТЕ 
However a volume much smaller than 0,6 cubic meters would be desir- 
able if accurate porosity predictions are to be made, It is pos- 
sible to limit the volume by insulating it from its surrounding 
environment (inside the corer barrel) or by focusing an electrode 
array (outside the corer barrel) in such a way as to direct the cur- 


rent patterns to a confined volume, 





IN SITU MEASURING SYSTEM 


The Electrical Resistivity Measuring System was built in the 
Ocean/Civil Engineering Laboratory of the University of Rhode Island 
for the DOSP, It consists of three components: electrode arrays, 
an electronic measuring circuit and a FM multiplex telemetry data 


link, 


Electrode Arrays 


The Electrical Resistivity Measuring System for the DOSP has 
two types of electrode arrays: the electrode array used inside the 
corer tube and the electrode array used outside the corer tube. The 
electrode array used inside the corer tube (see Figure 36) consists 
of four brass ring electrodes vertically spaced along the inner wall 
of the corer tube. The corer tube of the DOSP is a continuous five 
foot section of transparent acrylic tubing (0.D. 3 in., I.D. 2 3/4 in.). 
An ditional 14 inches of acrylic tubing (0.D. 3 1/4 in., I.D. 3 in.) 
is fused to the lower end of the corer tube. This allows the adapta- 
tion of a corer cutter-catcher assembly and provides additional struc- 
tural strength in the vicinity of the inner electrode array. The 
placement of the inner corer tube electrodes was difficult because, 
the inner corer tube electrodes had to be placed flush against the 
inner corer wall, Electrode extension into the sediment causes 
sample disturbance when the sediment sample is taken while recessing 


the electrodes reduces the electrical contact with the medium sampled. 
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Figure 36. Nose of the DOSP Corer Tube with the Inner Corer 
Ring Electrode Array and the Corer Cutter-Catcher 
Assembly. 
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In addition the inner ring electrodes require extra acrylic tubing 
material for structural support. The apparent resistivity (ra? of 
a marine sediment sample through the corer tube is calculated by 


the following formula: 


N Kal, 
where: 
d is the inner diameter of the corer tube; 
a, the separation between the inner two potential electrode 
rings; 
AV, the potential difference between the voltage electrodes; 
I che current; 


К, the geometric coefficient for the particular electrode 
2 `~ ` 
та 


array is equal to ER * 


Equipotential surfaces for a vertical cross section of the inner 
corer tube electrode array are shown for a homogeneous medium in 
Figure 37. From the equipotential surfaces an estimate of the volume 
contributing to 80 percent of the total resistance involved in a 
resistivity measurement is 580 cubic centimeters, 

The electrode array used outside the corer tube is based on 
the principle used in the oil well logging device called the 
Microlateralog. The Microlateralog was introduced by Dou of the 
Schlumberger Well Surveying Corporation, The device was designed to 


indicate the resistivity of a small volume of material immediately 


in front of the electrode array. The Microlateralog (see Figure 38) 
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Microlateralog Pad Electrode Array and the 80 
Percent Equipotential Line. (Shaded area is 
approximately 80 percent of the total ohmic 
drop of potential). 
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consists of a center current electrode and three circular ring 
electrodes. The largest or outer ring is the other current elec- 
trode and the inner two rings are the voltage electrodes. The 
electrodes are fabricated out of non-corrosive monel and embedded 
in an insulating supporting material, The small pad can be easily 
mounted either outside the corer barrel just above the corer 
cutter-catcher assembly, or on one of the DOSP probe legs. The 
vertical cross section showing the direction and shape of the cur- 
rent lines for a homogeneous medium are also shown in Figure 38. 
The dark equipotential surface represented in Figure 38 is such 
that 80 percent of the potential drop in the current pattern occurs 
between the electrode pad and this surface. Consequently, the 
volume of material confined within the pattern extends to a distance 
of approximately 3.5 inches from the electrode pad and contributes 
80 percent of the total resistance involved in the resistivity 
Measurement, This volume for a homogeneous medium is approximately 
380 cubic centimeters, 

Since both the inner and outer corer barrel electrode arrays 
require wire connection along the corer tube to the measuring instru- 
mentation, a flat, 9 conductor, polyester insulated, transmission 
cable, manufactured by the Burndy Tape Cable Div., 15 Linden Park, 
Rochester, N. Y., is to transmit the current from the signal 
conditioner to the current electrodes and to send the potential dif- 
ference from the voltage electrodes back to the signal conditioner 
(see Figure 39). Four of the conductors are signal cables and five 


are ground cables; each signal cable is separated by a ground cable, 
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Photograph of the Inner Corer Ring Electrode Array 


Figure 39, 


Cable Hook Up. 
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The flat transmission cable thickness is 0.0085 inches. The cable 
was fastened о the corer tube with three successive thin filmed 
coatings of ^ cement, With this configuration the connecting 
wiring to the electrode arrays does not hinder corer penetration 
nor will the wiring break due to abrasive materials, When the 
Microlateralog pad was placed on the corer tube, the same type of 
flat transmission cable was used for electrode hookup. 

In addition a miniature guard mode electrode array, described 
in the Literature Review section, was fabricated for laboratory 
testing. Seven thin platinum wire electrode rings were vertically 
spaced on a plastic probe five inches long and 1/4 inch in diameter, 
Thin wire connecting the platinum ring electrodes to the electronic 
circuit was embedded in the plastic probe head. The miniature guard 
mode electrode array probe was not tested at sea since the plastic 
material used in fabrication was not considered sturdy enough to 
withstand sea tests, 

A photograph of a Microlateralog pad and the plastic probe 
with a guard mode electrode array is shown in Figure 40. A photo- 
graph of the inner corer ring electrode array is shown in Figure 41, 
Any of the three electrode arrays can be used with the electronic 
circuit, 

The acrylic corer tube with its inner brass electrode rings 
was fabricated by Plastic Fabricators Inc., East Greenwich, R. 1. 
The Microlateralog pad and guard mode probe used to measure resisti- 
vity outside the corer barrel were fabricated at the Naval Underwater 


Systems Center, New London, Conn, 
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Photograph of the Microlateralog Pad Electrode 


Figure 40. 
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Figure 41. Photograph of the Inner Corer Ring 
Electrode Array. 
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Electronic Measuring Component 


The Electronic Measuring Component consists of an electrical 
conductivity signal conditioner, DC/AC inverter, a constant voltage 
transformer and a thermistor assembly. All of the above components 
are located in a pressure case.(see Figure 42), m the exception 
. of the thermistor assembly which is located in a separate housing 
that is threaded into the upper end-cap of the pressure case, 

The pressure case with the electronics will be mounted on the 
DOSP. The pressure case is cylindrical with an inside length of 
24 inches and inside diameter of 6 inches, The top end-cap is tapped 
to allow the threading of the electrical connectors (manufactured by 
Vector Cables Corporation, 5615 Lawndale, Houston, Texas) and the 
thermistor housing. Through these electrical connectors the elec- 
tronic кти гас receives 12 volts, 1 amp. D.C. power from the 
DOSP batteries, In addition, signals are transmitted and received 
from electrode arrays through these connectors, and a D.C, data sig- 
nal is sent to the FM telemetry system for transmission to the sur- 
face. Both pressure case end-caps are removable and the electronic 
components are mounted on trays which are connected to the top end- 
cap. The pressure case wall thickness is 1/2 inch for the cylinder 
and 1-1/2 inches for the end-caps. sine fabrication is ee 
6061-T6 aluminum, a thickness diameter ratio of 0.083 will permit a 
fydrostatic pressure of 5000 psig with a safety factor of 1.30. The 
end-caps have the same pressure capabilities. The pressure case is 
tested repeatedly to 3000 psig in the University of Rhode Island 


hydrostatic test tank before actual sea tests. The procedure for 
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Pressure Case and Electronic Measuring Components: 


Figure 42. 


Signal Conditioner, DC/AC Inverter, Constant Voltage 


Transformer and a Thermistor Assembly. 
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safety testing pressure cases on the DOSP is more thoroughly 
discussed in Lewis.” 

The conductivity signal conditioner, manufactured by 
Honeywell Inc.; Fort Washington, Pa., (Model description number 
552022-1002-10Q-000), is a solid state null balance A.C. amplifier. 
It supplies the voltage to drive the sensor array, and compares the 
resultant output of the sensor voltage electrodes with a stable ref- 


erence voltage 


s А differential voltage is used to regulate the 

? 

drive current. The signal conditioner produces an output linearly 
proportional to the conductivity of the measured medium. The out- 
put signal range Т O to 5 volts D.C. and is sent to the FM tele- 
metry system. 

The signal conditioner operates as a potentiometer. Ап А.С, 
current I. is applied through two sensor current electrodes, Cl and 
C2 (see Figure 43). The resulting voltage drop, Vos across two 
measuring electrodes, Ml and M2, is compared with a fixed A.C. ref- 
erence voltage, Ув. When they rcm the resulting error voltage, 
Vos alters the sensor drive current I.» The amplifier system then 
changes its output to maintain уш equal to Vee Since ү, is held 
constant, the drive current (т) is (by Ohm's Law) inversely pro- 
| portional to the medium's resistance and directly proportional to 
its reciprocal conductivity. 

Because of the large measuring error which would be intro- 
duced by polarization of the sensor electrodes, A.C. rather than D.C. 


current is employed. The voltage required for the signal conditioner 


is from 107 to 127 volts A.C. at 50 or 60 Hz. This is the same 
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Thermistor 
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Figure 43. Skematic Diagram of Electronic Measuring System 
and the Inner Corer Ring Electrode Array. 
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frequency supplied to the sensor electrodes and is high enough to 
prevent polarization. The A.C. power requirement of the signal 
conditioner makes an inverter necessary. The inverter, a 12 volt 
D.C. input,square wave 60 Hz output, is manufactured by Fork 
Standards, Іпс., West Chicago, Ill. (Model number 32A057d). Ве- 
fore the signal is sent on to the signal conditioner, a constant 
voltage transformer (produced by Solar Electric Co., Chicago, 111, 
catalog number 30498), is used to smoothe the signal out. The 
inverter and constant voltage transformer allow the use of the 

` DOSP battery power supply. This permits the Electrical Resistivity 
Measuring System to be used on the ocean floor in water depths of 
5000 feet, 

The thermistor assembly and housing (type number H66), is 
manufactured by Fenwal Electronics Inc., Framingham, Mass. The 
nominal pressure rating of the housing is 10,000 psig. The ther- 
mistor automatically corrects the conductivity readings to a reading 
for 25° Centigrade. The thermistor is connected to an adjustable 
compensation network in the signal ce which alters the ref- 
erence voltage to match sensor output versus input voltage charge 
with the medium's temperature. The range of adjustment in the sig- 
nal conditioner is from 1,4 to 2.5 percent per degree Centigrade. 
For most marine conditions, 2.5 percent per degree Centigrade is 
chosen since this is the temperature coefficient usually used for 
sea water, For shallow water tests the thermistor was embedded into 
fhe acrylic corer barrel next to the brass electrodes, The ther- 


mistor was wired to the cable wire and sealed with epoxy cement. 
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FM Multiplex Telemetry Component 


The development of a telemetry component for the DOSP was 
conducted by Barry i under contract with the Naval Underwater Systems 
Center, The system was designed to command the DOSP from a remote 
location on a surface ship, and receive and transmit data signals 
from the DOSP in ocean depths to 6000 feet. All signals are sent 
and received via coaxial cable by means of a frequency multiplexed 
(FM) technique. 

The FM telemetry system for the DOSP can be divided into three 
distinct spectra: touch tone signals (25 to 115 Hz), the slowly 
varying channels (370 to 15.5 kHz), and the acoustic channels 
(150 to 450 kHz). Resistivity is one of the slowly varying data 
channels with a channel center frequency of 1300 Hz. 

The Electrical Resistivity Measuring System's signal condi- 
tioner provides the necessary voltage level (0 to 5 volts) to drive 
the voltage controlled oscillator of the telemetry system, After 
the signals have been frequency modulated they are multiplexed with 
acoustic data and touch tone commands. This technique allows all 


data to be transmitted simultaneously. 





PROCEDURE 


Laboratory Tests 


Laboratory tests were conducted after the Electrical Resis- 
tivity Measuring System was fabricated in order to test its perform- 
ance prior to sea tests. For analysis and calibration all three | 
components (electrode array, electronic circuits and FM telemetry 
data link) were connected together and data recorded, The full 
range тыг (0 со 5 Mits D.C.) signal that is accepted by the 
FM telemetry data link was supplied by the signal conditioner, This 
full range was obtained by increasing the salinity of the solution 
that the electrode array was sampling, The full scale conductivity 
range that may be sampled is from 0 to 60,000 y mho/cm. 

The inner corer ring electrode array was ee by 
filling the corer tube with standard solutions of known conductivity 
while the Microlateralog pad and miniature guard mode electrode ar- 
rays were immersed in the standard solutions, The voltage output 
for each electrode array in the standard solutions was read on the 
Signal conditioner, 

When the Electrical Resistivity Measuring System was used 
independently of the DOSP for testing in the laboratory or in 
shallow water, the FM multiplex telemetry component was no longer 
required. Under these conditions power was supplied either by D.C. 


battery or by standard 115 volts, A.C. 60 Hz that have passed the 
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inverter. Thus wiring to the electrode array and a volt meter to 
read the signal conditioner output voltage were the only items 
required to make the conductivity measurements, 

Ottawa sand and Providence silt were used in the laboratory 
tests in order to compare the results of measurements obtained with 
the different electrode arrays with the Т from earlier laboratory 
experiments, The procedure for testing the inner corer ring elec 
trode array was as follows: a sodium chloride solution of known 
conductivity was first poured into the corer tube which had been 
capped to prevent leakage. A measurement of the conductivity of 
this solution (к) was then taken in the corer tube. A known weight 
of Ottawa sand was then added slowly. The sand settled out rapidly 
in a loosely packed configuration, The length and the conductivity 
of the sand-water matrix (к) were measured and recorded for later 
porosity calculations. On completion of this measurement, the corer 
tube was tapped five times to densify the sand. Conductivity meas- 
urements were again taken and the length of the sand-water matrix 
was again cd: This procedure was repeated until the conducti- 
vity and length measurements became constant, The sand was then in 
a closely packed configuration, 

The Microlateralog (pad) and the ie guard mode elec- | 
trode array (probe) was tested by using a large acrylic tube chamber 
five feet long and one foot in diameter filled with 40 inches of 
Providence silt and 10 inches of saline water (30 parts per thousand). 
The pad and the probe were placed, one at a Si in the water direct- 


ly above the sediment interface and the water conductivity (к) was 
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measured. In turn both the pad and probe were lowered into the silt 
and at six inch increments the conductivity of the sediment (к) was 
measured. The, larger size of the pad required greater force to 
penetrate the нее; consequently it was only inserted to a depth 
of one foot. The probe was inserted the entire length of the silt 
column, Laboratory analysis of the marine sediments was conducted 
in order to compare the TD with those of the in situ sea tests, 
Water content, void ratio, shear strength, specific gravity, organic 
content, pH, grain size analysis, liquid limit and plastic limit 
were determined. In addition X-ray diffractograms were obtained for 
the clay-silt size fraction of each marine sample. Details of the 
above analysis are presented in Appendix C, 

The interstitial water salinity was obtained after the water 
content of the marine sample was analyzed. The dried sample was 
mixed with twice as much distilled water as was originally present, 
The slurry was repeatedly mixed and allowed to settle out, then a 
conductivity measurement of the solution was taken with a Beckman 
conductivity meter. Since the number of ions in solution remained 
constant and conductivity changes are linear with dilution change, 
the original conductivity of the interstitial water can be determined 
by doubling the conductivity reading of the dilute solution. These 
measurements were compared with in situ and laboratory measurements 
of bottom water salinity. 

Part of the marine sediment was diluted into a slurry using 
bottom water obtained from the core sites and then was re-consolidated 


to pressures of 0.25, 1.60, 4.00 and 16.00 psig in the redesigned 
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laboratory equipment described earlier. This allowed a Formation 
Factor-porosity канен нер to be obtained. Comparisons were also 
made with both the in situ Formation Factor and porosity values 
obtained from RT analysis of the marine sediment from the 


Electrical Resistivity Measuring System corer, 


Sea Tests 
Shallow water sea tests were conducted in Narragansett Bay 
at the following locations, which were chosen for convenience and 
for the bottom type variety found there; 
о : O : 
Station A- 41 34.1 N - 71 25.4 М 
O : | O : 
Station B- 41 29.5 N - 71 24.7 W 
о ' O 5 
Station C - 41 22.78 - 71 30.6 М 
The inner corer ring electrode array was tested at Station A. 
Station A had a sandy silt bottom with a 7 to 18 percent clay frac- 
tion, Since the depth of the water was approximately 20 feet, the 
corer barrel was attached to two ten foot sections of 14 inch dia- 
meter galvanized pipe. The procedure at this station was as follows: 
The research vessel was moored with bow and stern anchors and the 
Station location was confirmed by a visual fix. The salinity and 
temperature of the water column were measured by a Beckman model 
RS5-3 in situ salinometer-temperature meter. While the corer barrel 
was lowered slowly to the bottom the conductivity was continuously 
measured. Since the FM telemetry component was not needed in shal- 
low water, a voltmeter was used to read voltage output directly from 


the signal conditioner. As the corer barrel penetrated the bottom, 
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conductivity was measured every six inches. The corer barrel was 
then pulled out and contínuous conductivity readings were made un- 
til the core was stored on deck. Conductivity readings were again 
taken in the laboratory. 

Station B, the site of the future Department of Ocean Engineer- 
ing underwater laboratory, had a silty em bottom with a 
4 to 9 percent clay fraction. The depth of the water was approxi- 
mately 40 feet. Both the inner corer ring electrode array and the 
Microlateralog pad array were tested at this station. A six foot 
section of one inch diameter galvanized pipe was securely mounted to 
the Microlateralog pad which was then carried down, positioned, and 
inserted into the bottom by two trained scuba divers, The ring 
electrode array in the corer barrel was also taken down by the scuba 
divers. An attempt was made to systematically survey the area, 
Using the pattern shown in Figure 44, Microlateralog penetrations 
were made at radial distances of 25 and 50 feet from where the corer 
was inserted. The numbers circled in Figure 44 were the positions 
that were measured before the scuba divers' air supply was depleted. 
The procedure was the same as that used for Station A, except that 
the scuba divers handled the corer barrel and the Microlateralog pad. 
White tape was placed at six inch intervals on the corer barrel (see 
Figure 45) and every 12 inches on the Microlateralog pad to allow 
visual observation of penetration depth. The positions had been pre- 
marked with weighted cans. The Е. held the pad in the water 
prior to penetration to enable a bottom water conductivity measure- 


ment. The Microlateralog pad was inserted into the bottom at 
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CORE TAKEN 
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Figure 44. Survey Pattern Planned for Narragansett Bay 
Station B. 
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Corer Barrel with White Tape Marks 





Figure 45. 
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one foot increments. A communication line notified the divers (by 

a number of pulls) that the measurement was complete, and then the 
probe was inserted an additional foot. Since the Microlateralog pad 
did not have a thermistor, it was assumed that the temperature of 
the sediment was the same АЕ the temperature of the bottom water, 
From the later temperature measurement, the conductivity readings 
were corrected to 25° Centigrade, 

Station C, a shallow water area near Succotash Point, 
Jerusalem, Ra ES a sandy bottom, The water was 0 to 5 feet 
deep. Both the inner corer ring electrode array and the Micro- 
lateralog pad array were tested at Station C. The pad E. attached 
to the six foot galvanized pipe used at Station B. The procedure 
used was the same as for the other two stations, except that the 
'instrumentation was left on shore, the corer barrel and the pad were 
carried out by a wader, and measurements were made at random, The 
sand bottom was very hard, consequently the core obtained was ap- 
proximately one foot and the Microlateralog pad could only be in- 


serted to 6 and 12 inches. 





RESULTS AND DISCUSSION 


Laboratory Tests 


The results of the calibration of the total system (all three 
components) are presented in Figure 46, Figure 46 shows that the 
0 to 5 volt D.C. output signal which was sent through the FM tele- 
metry component and 6000 feet of coaxial cable was transformed line- 
arly to a range of -1.5 to 1.5 volts D.C. Calibration curves for 
the inner corer ring electrode array, the Microlateralog pad array 
and miniature guard mode array are shown in Figure 47, These curves 
indicate that, for each electrode array, an individual linear rela- 
tionship exists between the conductivity of the medium sampled and 
the signal conditioner output voltage. Consequently ЖА ЕГП 7 сап 


be written as: 


Im Ke um 
18 о ү | v 
с m m | 


where D is the signal conditioner output voltage when the electrode E 
array is sampling the continuous phase (water), 

and Ya is the signal conditioner output voltage when the electrode 
array is sampling the mixture (sediment). 

Therefore if these two voltages represented in Equation 7 are meas- 

ured with temperature compensation, the Formation Factor of a sedi- 


ment can be readily obtained, 


168 





Recorded Level (volts) 


5 а 5 


VCO input (volts) 
or ` 
DC voltage output 
5 of signal conditioner 


Figure 46. Calibration of FM Telemetry Channel showing 
Voltage Input from the Signal Conditioner and 
the Response Recorded Through FM Telemetry 
Channel. 
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The results of laboratory testing of the inner corer ring 
electrode array with Ottawa sand compared well with the results 
of the laboratory equipment described earlier (see Figure 48). The 
data from the inner corer ring electrode array are shown to be with- 
in the + 2 percent porosity error lines of the line of best fit 
determined earlier for Ottawa sand by 47 Formation Factor-porosity 
tests. 

In Figure 49 Formation Factors obtained with the miniature 
guard mode electrode array and the predicted porosity (obtained 
from Figure 22) are plotted versus depth for Providence silt. A 
sediment sample was obtained from the surface, the center, and the 
bottom of the Providence silt column for porosity determination. As 
shown in Figure 49 the value agrees well with the predicted уай е. 
Three measurements were taken with the Microlateralog pad array 
probe to a depth of one foot. These three Formation Factors along 
with the predicted porosity are also shown in Figure 49 as a function: 


of depth, and agree well with other data on Providence silt, 


Sea Tests 

In Figure 50 ten in situ Formation Factors for Station A and 
five in situ Formation Factors for Station B are К... as a func- 
tion of the porosity determined in the laboratory. The lines of 
best fit which were determined from 28 measurements for Station A 
and ll measurements for Station B by the redesigned laboratory equip- 
ment are also shown in this figure. Two thirds of the in situ Forma- 


tion Factors fell within + 2 percent of the lines of best fit while 


FORMATION FACTOR 
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Figure 48. Formation Factor-Porosity Data Points of Ottawa 


Sand Obtained with the Inner Corer Ring Electrode 
Array and the Line of Best Fit from Earlier Laboratory 
Data of Ottawa Sand. 
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X- ACTUAL POROSITY MEASUREMENT 

o- FORMATION FACTOR (GUARD MODE) 
E)-FORMATION FACTOR (MICROLATERALOG) 
A-PREDICTED POROSITY 

V-PREDICTED POROSITY 


Formation Factors Versus Depth Obtained with the 
Miniature Guard Mode Electrode Array, the Corres- 
ponding Predicted Porosity and the Measured Values 
of Porosity. | 


173 


174 


2% POROSITY 
ERROR LINE 


FORMATION FACTOR 





A-A LINE OF BEST FIT (STATION A) 
OBTAINED FROM REDESIGNED Equipment A 
©-CORER | STATION A DATA 
[J-CORER 2 STATION A DATA 


B-B LINE OF BEST FIT (STATION B) 
OBTAINED FROM REDESIGNED EQUIPMENT 
A-CORER | STATION B DATA ` 





30 40 50 60 70 80 90 00 
POROSITY (%) 


Figure 50. In-situ Formation Factors Obtained with the Inner 
Corer Ring Electrode Array Versus the Measured 
Laboratory Porosity Values for Station A and B 
and the Lines of Best Fit for Stations A and B 
Obtained From the’ Redesigned Laboratory Equipment. 
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one third were within + 4 percent, This same information is 
presented versus depth of the core in Figure 51, 

There are several factors that could have caused this vari- 
ation. The variation was not caused by assuming the interstitial 
water conductivity was the same as the bottom water conductivity, 
since as shown in Figure 52 the interstitial water salinity deter- . 
mined in the laboratory on the cores taken compare well with the 
in situ bottom er salinities. For example the interstitial 
water salinity of the second core taken from Station A ranged from 
31.0 to 31.8 parts per thousand while the values obtained by the 
inner corer ring electrode array, the in situ salinometer and the 
laboratory conductivity bridge were 31.2, 31.1, and 31.0 respective- 
ly. These results agree with the work of Ruiter bees = and 
К..." and justify the assumption used by Kermabon et al” that 
Formation Factors could be determined by using bottom water salinity 
in place of interstitial water salinity. 

The fact that one third of the Formation Factors fell outside 
the + 2 percent error lines would not appear to be caused by the 
equipment since the laboratory tests with this equipment did not 
have such a wide variation, The presence of shells is another pos- 
sible cause of variation, ГАНЕ sample disturbance before the 
porosity of the sample was determined in the laboratory appears to 
be the most reasonable explanation. Although disturbance was intended 
to be monitored by continuous conductivity readings during core re- 
trieval, this did not occur at Stations A and B due to an error. 


introduced by the temperature thermistor. The thermistor was de- 
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signed to automatically correct the conductivity readings to 

25° Centigrade, and the error was introduced whenever the temper- 
ature sensing element embedded. in the acrylic corer barrel was dif- 
ferent from the temperature of the sediment. This was observed as 
the cores were retrieved from Stations A and B after the cores were 
removed from the refrigerator in the laboratory, This error could. 
be eliminated and disturbance adequately monitored in the following 
wayss | 

1. Measure the conductivity of the core sample after it has 
reached room temperature, 

2. Measure the conductivity of the core sample in the refrig- 
erator (or ship's reefer) after the sediment has reached 
the refrigerator temperature, 

3. Bypass the thermistor in the electronic circuit; measure 
the temperature of the sample and correct the conductivi- 
ty measurement to 25 degrees Miitigrade. 

If the thermistor was bypassed and the temperature corrected 
to 25 degrees Centigrade by taking the temperature manually, the 
conductivity changes upon retrieval could be considered due to 
ШЕ ace. If this were repeated in the laboratory, any changes 
in the conductivity readings would be most likely due to disturbance. 
Disturbance indicated by a settling of about one inch was observed 
on the core from Station В, 

Proof that disturbance of the samples before the porosity was 
termined could have caused the wider variation in porosity predic- 


tion of the sea tests is shown in Figure 53, The sediment column 
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Figure 53. Formation Factor-Porosity Relationship at Different 
Stages Of Sediment Disturbance in the Inner Corer Ring 
Electrode Array and the Line of Best Fit of the Same 
Sediment Obtained by Using the Redesigned Laboratory 
Equipment. 
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height decreased 2,2 cm. and the porosity decreased from 41 percent 
to 38 percent after the corer from Station C was hit and dropped. 
As seen ín Figure 53, Formation Factors that are not corrected а 
this sample disturbance lie on the + 4 percent error line. The line 
of best fit was determined for 13 measurements by the redesigned 
laboratory equipment. Five Formation Factor-porosity values were 
determined at different stages of disturbance after the core stood 
for 24 hours at room temperature, Since the sediment column height 
was taken at each disturbance stage, the porosity was determined by 
using the diameter of the corer barrel and the weight of the soil. 
A — calculation is presented in Appendix C, Figure 53 shows 
that if this INES is taken into account the porosity predic- 
tions lie within + 2 percent of the line of us fit. This shows 
that the inner corer ring electrode detected the disturbance when 
the error due to temperature differences was eliminated. When the 
porosity determinations were correlated with resistivity measure- 
ments taken after disturbance, the Formation Factors agreed well with 
the line of best fit. 

An important advantage of the Electrical Resistivity Measur- 
ing System corer was discovered when the core cutter-catcher assembly 
did not work and the core was lost at Station B. The loss of the 
core was immediately detected by the fact that the conductivity read- 
ing changed to that of the surrounding sea water, Thus in deep sea 
operations another core can be taken without bringing the corer to 
the surface, 


The inner corex ring electrode array can be used to ^;  vmine 


| 
| 
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the salinity of the water column. At Station A the conductivity 
readings of the water column indicated a salinity of 30.5 parts per 
thousand while the in situ salinometer read 30.2 parts per thousand, 

The predicted porosity determined by the Formation Factor 
measurements made by the Microlateralog pad was obtained by using 
the line of best fit for Station B shown in Figure 50. These pre- : 
dicted porosity values are shown as a function of depth (to 5 feet) 
for six positions at Station B in Figure 54. For comparison, poros- 
ity as determined in the laboratory, of the two foot core obtained 
from the: Electrical Resistivity Measuring System corer is also shown 
in Figure 54. The porosity Pues do not agree as well as expected, 
For example at a depth of one foot, the predicted porosity value of 
position one (25 feet from the core) was 72 percent while the Jaa 
determined on the core was 66 percent, However the predicted poros- 
ity determined by the Formation Factor measurements made by the 
Microlateralog pad for six positions at Station C and the line of 
best fit for Station C (see Figure 53) was in good agreement with the 
porosity determined in the laboratory. The average Tree value 
was 42 percent for six inch pad penetration and 40 percent for one foot 
penetration while the porosity values determined in the laboratory 
ranged from 38 to 41 percent, Since the penetration was limited to 
one foot in the hard sandy bottom of Station C, a vibrating or pneu- 
matic hammer similar to Borowczyk's and Krolikowski's” Universal 
Probe would have to be used for deeper penetrations, 

The Formation Factors determined in situ and the laboratory 


determinations of porosity, water content, void ratio, cone penetra- 
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meter, organic content, Atterberg limits and pH data on the two 
cores from Station A and the one core from Station B are shown in 
Figures 55 and 56 respectively. The in situ Formation Factor varied 
inversely with the porosity К water contents. These results agree 
with the findings on the marine sediments discussed in the first 
section of this paper. Although the laboratory analysis presented 
in these figures help define the sediment type at these stations, no 
other soil property correlations can be made from this data. If a 
larger number of cores were taken, it might be possible to correlate 
resistivity with another soil property such as shear strength, bulk 


density and acoustic properties. 
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Figure 55. In Situ Formation Factors Porosity, Water Content, Void 
Ratio, Cone Penetrameter, Organic Content, Atterberg 
Limits, and pH of Two Cores from Station A. (Organic 
content, Atterberg Limits and pil were determined from 
top, bottom and middle core sections) 
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CONCLUSIONS 


1. The Electrical Resistivity Measuring System that was designed 
Яна fabricated was successfully tested in the laboratory and in 
Narragansett Bay. The system was tested by connecting the electrode 
array, the electronic circuits, and the FM telemetry data link with 
6000 feet of coaxial cable in order to simulate the requirements for 
deep ocean testing on the Deep Ocean Sediment Probe (DOSP). Three 
electrode arrays were fabricated and tested: the inner corer ring 
electrode array, the Microlateralog pad electrode array and the mini- 
ature guard mode electrode array. The predicted porosity values ob- 
tained with the system in the laboratory were in good agreement 

(+ 2 percent of the lines of best fit) with values determined earlier 
with the redesigned laboratory equipment. The system was tested in 
the shallow water of three locations in Narragansett Bay. Using the 
inner corer ring electrode array most of the data obtained was in 
good agreement with the data u later in the laboratory with 
the redesigned laboratory equipment. The data that was not in good 
agreement (+ 4 percent of the line of best fit) could have been due 
to disturbance of the sample before the porosity was determined in 
the laboratory. The Microlateralog pad electrode array needs 
further testing. Though the porosity prediction using the line of 
best fit was in good agreement with laboratory porosity values at 


Station C, a wider variation was found at Station B, 
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2. The conductivity readings of the inner corer ring electrode 
array on the water column and the bottom sea water converted to a 
salinity very similar to the readings ae an in situ salinometer, 

30.5 and 30.2 parts per thousand respectively, The salinity of the 
bottom water as determined by the inner corer ring electrode array 
was in good agreement with the interstitial water salinity deter- 
mined in the laboratory, 31.2 and 31.1 parts per thousand respec- 
tively. 

3. The Electrical Resistivity Measuring System was found to be 

an excellent monitoring system, The System indicated that the bottom 
was reached by showing a marked change in conductivity readings. In 
the same way the loss of a core was immediately detected, This abil- 
ity would be especially valuable for a deep ocean corer. Tests in 
the laboratory with a sediment in the Electrical Rewer ity Meas- 
uring System corer showed that the resistivity measurements detected 
sample disturbance, 

4. The Electrical Resistivity Measuring System costs approximately 
$1,500, not including the FM telemetry data link component which is 
shared by other systems aboard the DOSP, This relatively low cost 
does not include the 50 man hours required to connect the system, but 
does include all the electronic circuits, electrode arrays (including 
the corer barrel), pressure housing, thermistor housing, electrical 
connectors and wiring, Since the electronics of the Electrical Resis- 
tivity Measuring System and the redesigned laboratory equipment are 
the same, module components of these systems can be interchanged. 


This duplicity affords quick back-up spare parts whenever necessary. 
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The Electrical КОО System and its components 

could be adapted to a probe, coring device or any deep sea plat- 
form. It could be used ín the deep ocean or in shallow coastal 
waters, If hard sandy bottom is encountered, a vibrating probe 


would be required for deep penetration, 





FURTHER STUDIES 


1. The results of this study indicate that electrical resisti- 
vity measurements could be used to determine the cation exchange 
capacity of individual clays. For example clays with a high cation 
exchange capacity could be expected to have Formation Factors of less 
than one at low interstitial water conductivities. Clays with a low 
cation exchange capacity would always have Formation Factors greater 
than one and the Formation Factors would remain se regardless 
of the interstitial water conductivity (assuming porosity remained 
the same). The Formation Factors of the individual clays could be 
correlated with the cation exchange capacity obtained by standard 
methods. Once established, the electrical resistivity method would 
be easier and less time consuming than the standard method, 

2. The development of an induction type Electrical Resistivity 
Measuring System should be pursued for in situ sediment measurements, 
This type of measuring system would eliminate direct metal electrode 
contact of the sediment and could result in greater accuracy. 

3. The combination of the Resistivity, Acoustic, Nuclear, Thermal 
and Geochemical measuring systems on the DOSP offer a unique opportu- 
nity to further study and determine accurately in situ properties 
such as porosity, water content, sound velocity and bulk wet density, 
In addition the ability of the Electrical Resistivity Measuring Sys- 


tem to monitor such situations as touch down, core loss, and sample 
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disturbance make the system even more valuable for deep water. 

These side assets should be further studied in future deep sea 

tests. 

4, The use of the Electrical Resistivity Measuring System in 
pollution detection should be studied further. For example a 
continuously monitored water pollution detection device could be 
simply attached to a buoy. The cell would consist of two small 

tubes each having an inner ring electrode array. One tube contain- 
ing water of desired qualities would be sealed, The бекет tube, 
adjacent' to the first, would be a flow-through tube sensing the 
environmental water, The instrumentation would work on the Formation 
Factor concept, i.e. if the environmental water exactly met standards 
the Formation Factor would be one. Any large variation from unity 
would give warning of a possible polluting situation, 

Perhaps a similar Formation Factor device could be used to 
measure the cleansing or infusion of pollutants in bottom sediments. 
These residual pollution effects in bottom sediments could be moni- 
tored continuously and over long periods of time, 

5, New electrode arrays similar to the probe described by 

Bouma et gu? could be tested for a wide range of uses. For example, 
these probes could be inserted to considerable depths (20 feet) at 
various positions in an estuary such as Narragansett Bay, and 

seasonal variations in porosity could be detected, 

This same electrode array placed on the legs of an oil rig 
platform or Texas tower could serve as a warning system for such 


phenomena as bottom erosion and scour. 


192. 


6. The ability of electrical resistivity measurements to predict 
other submarine soil properties in situ should be pursued. For ex- 
ample, undrained shear strength measured in situ by vane shear or a 
penetrameter could be correlated with electrical resistivity measure- 
ments, Also relationships could be developed between electrical re- 
sistivity measurements and acoustic properties of soils or geochemical 
properties such as Eh and pH. 

#% The possibility of attaching an electrical resistivity elec- 
trode array to a telescoping or vibrating deep ocean floor penetra- 
ting device (50 feet or more) merits further study since properties 
such as porosity, water content and density could be predicted to 
these depths, This knowledge would be very useful to ocean engineer- 
ing projects such as ocean floor foundations. 

8. A sample could be monitored by the Electrical Resistivity 
Measuring System while it is returned to the surface from the deep 
ocean if a second temperature sensing element was provided. Since a 
thermistor is located in the tip of one of the DOSP's probes, a sec- 
ond temperature of the water column could be used to adjust the con- 
ductivity measurement to a correct value, and the disturbance of a 
sample during retrieval could be detected. This method of monitoring 
Afeturbance could be checked by the nuclear gauge suggested by 

Rose and aT or an acoustic dor | 

9, Electrical resistivity measurements could be used to determine 
sediment grain size, This could be accomplished by fabricating a 

30 inch long, one inch diameter plastic tube with three or four sets 


of inner ring electrode arrays placed vertically along the inner wall 
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of the tube, 

A slurry containing a known amount of soil, water and de- 
flocculating agent would be added to the cell tube and spec Canes 
resistivity measurements would be made at definite time increments, 
As the soil ата settled, the resistivity of the lower array 

would increase while the resistivity of the upper arrays would de- 
crease. If the same assumptions used in a hydrometer analysis are 
applied to this method, the electrical resistivity measurement could 
then be related to particle size. This method would be easier than 
the hydrometer method since the measurements could be made. electri- 
cally and a graphical record could bs made automatically. 
10. Electrical resistivity measurements could be used to monitor 
consolidation in situ. This could be accomplished by placing an 
electrode array under an object or structure with the electrodes’ 
current field in the direction of the sediment. Changes in ne 
porosity due to the weight of the structure would be detected by the 
electrical resistivity measurements, Since the cost of electrode 
arrays is low, several could be placed to monitor changes both 


continuously and periodically. 
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APPENDIX A 


DETAILED DEVELOPMENT OF SOME EQUATIONS 


1. Development of “Parallel” Equation 

Assume that the medium consists of a series of plates, all par- 
allel to the direction of er flow, Then the plates may be thought 
of as individual resistors in parallel. (Some of the plates are sheets 
of interstitial water.) To find the total conductivity, the individual 


conductivities are summed, but weighted by their volumes: 


Kam = EeVe + Kala / 
where: 
Va 18 the volume of the mixture (sediment), 
V, is the volume of the continuous phase (water), 
Үз із the volume of the dispersed phase (soil). 


dividing by Va : 
E, = коп + Ka (1-п), 


where n is the porosity, and dividing by Ko ; 


Ku ge - 
=s kant K (l-n), 


or the Formation Factor (FF) is, 
re а 
FF = к. =n + Re (1-n). 
2, Development of "Series" Equation 


Assume medium consists of series of plates, all perpendicular to 


Al 


direction of current flow. To find the conductivity of the mixture, the 
individual resistivities are weighted according to their volumes and 


summed: 
ro r BUE тауа F 


dividing by Vn ; 


een т (1-n), 
dividing by r, , 
r Т 
= = д + = (len), 
c C 
Dr, 
FF = 


Қ, = n + R4 (SNE. 


3. Detailed Derivation of the Resistivity Equations for the Wenner and 
the Schlumberger Configurations of Electrodes 

The most direct approach to measurement of the value of sediment 
aggregate resistivity is by the four-electrode method. The .four-electrode 
system of measuring earth resistivity is quite old, having been developed 


by Wenner?” 


of the National Bureau of Standards in 1915. Referring to 
Figure A-1 a known current is passed between the outer two electrodes, 

and the difference in potential between the inner two electrodes is 
measured. The resistivity can be calculated from knowledge of the current, 
the potential difference, and the geometric arrangement of the electrodes, 


Assuming a point current source in а hamogenous medium, the Laplace 


equation may be written as: 


A2 
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If there is no angular dependence (i.e. the resistivity of the sediment 
is not, in general, perfectly isotropic) the Laplace equation can be 


written as 


( a“v/ar? ) + (2/R) (avV/eR) = 0 , 


/ 


where V is the electric potential at a distance R from the point source. 
The solution to this equation is well known and is 

У = ( С, /R) + С, : 
To determine the constants С, and C, , assume that as R extends to infin- 


ity, V tends to zero. Thus, C, = 0, and so 

у= с/в. | (1) 
By definition, the resistivity of a medium is 

r = -(A/I) (dV/dR) , (2) 


where r is the resistivity, I is the current, and A is the area through 


which it flows. Now from (1), 
dV/dR = CR (3) 


If it is assumed that a point is the source, then the current flows 
through a sphere surrounding the source, and the area of this sphere is 
simply ATR? . Substituting this value and equation (3) into equation 


(2), it becomes: 


A4 | 





r = (атт JIN (-с/в^) zw 4m C/I , 
or 
С = ак. 
Substitution of this value into (1) gives 
V = тв. (4) 


À current source implies a current sink, and the potential of such a 
sink is just the negative of (4). By superposition, the potential at a 


point p, V(p), is 
+ S 
Vip) = V (p) + V (p), 


+ z 
where V (p) and V (p) are the potentials at point p due to the current 


source and the current sink, respectively. Thus, 


Vip) = (Gi/4m) Gp RE), 


where Ry and Ro are the distances from the point p to the current source 
and sink, respectively. 

In practical situations, one measures the difference in potential 
between two points a finite distance from the source and the sink, as in 


the figure below, 


AS 


source point 1 point 2 sink 
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PTA 
ТЕЛУ 15 defined аз: 


AV = (potential at point 1) - (potential at point 2) . 


Then 


I Den l =- l = 1 
ДУ =з „ai nn ter, — nn e (5) 


The choice of the electrode spacing is somewhat arbitrary. The two most : 


common are the “Wenner” and the “Schlumberger” configurations. These are 


defined by: 
Wenner Schlumberger 
Ri a a 
К. 2а За 
К. да За 
R 4 a a 


where a is a convenient unit of spacing. 


Аб 





Equation (5) then reduces to 





Wenner г = 47ra(AV/I) 


Schlumberger г = 3fTa(AV/I) . 
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APPENDIX B 


DETAILS OF EQUIPMENT AND CALIBRATION 


The basic features of the laboratory and in situ equipment, with 
descriptive photographs, have been discussed in the text. Also calibra- 
tion procedures were discussed and in some cases calibration curves were 
presented. This section will present some detailed drawings of this 
equipment, block diagrams of electronics, some flow patterns developed 
by the electrode arrays and calibration curves not previously given in 


the text. : 


В, Drawings of Cells and Electrode Arrays 
A. Simple Cell Apparatus. 

Figures 7, 8 and 9 in the text are photographs of the simple cell 
apparatus. Figure B-1 is a drawing of the simple cell with two types of 
top electrodes used in testing. 

B. Modified Apparatus. 

Figures ll and 12 in the text are photographs of one sediment cell 
and chamber respectively. A detailed drawing of one cell tube is given 
in Figure B-2, Other components of the modified apparatus such as sum 
port frame, piston rod, piston head and bottom cap can be obtained from 
plan no. G-120 on file in the University of Rhode Island Engineering 
Instrument Shop. 

С. Redesigned Apparatus, 

A photograph of the redesigned apparatus is presented in Figure 

14, A drawing of the redesigned cell and pressures vessel is shown in 


Figure B-3, 
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Figure B-1 Simple Apparatus Test Cell with Variable Sample Length 
Electrode and Normal Top Electrode, 
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Figure B-2 Modified Apparatus Test Cell. 
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Figure B-3 Redesigned Apparatus Pressure Vessel and Test Cell. 
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В. Inner Corer Ring Electrode Array. 

In the text Figure 36 shows a descriptive о: the inner 
corer ring electrode array and Figure 39 is a photograph of the same 
array. A drawing giving the dimensions of this array is shown in Figure 
В-4. 

E. Microlateralog Electrode Array and the Miniature Guard Mode Elec- 
trode Array. 

Descriptive and dimensional detail is shown in the кинин. of 


these two arrays in Figure 40, , 


2. Block Diagrams of the Electronic Вано | 

А. Те EL S equipment used in the simple and modified apparatus 
can be seen in Figure 13. Figure B-5 gives the block diagram of this 
electronic equipment. 

B. The block diagram of the Honeywell Signal unes (model no. 
552022-1002-100-000) is shown in Figure B-6. The signal conditioner was 
the basic electronic component used in the redesigned laboratory apparatus 
and the in situ Electrical Resistivity Measuring System. 

С. Figure 42 in the text shows the electronic components of the Elec- 
trical Resistivity Measuring System. Figure B-7 is a block diagram of the 


system’s electronic components. 


3. Flow Patterns of Electrode Arrays. 

A. Simple cell electrode array for the current flow between top and 
bottom electrodes resulted in the flow pattern shown in Figure B-8. If 
the current was allowed to flow between the top and bottom platinum point 


electrodes the flow pattern shown in Figure B-9 developed. 
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Figure B-4 
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Figure B-5 Electronic Block Diagram of the Simple and Modified 
Apparatus. 


B7. 


UA e o dne 


313839408 


— — — — 


битом орча ys ru. 


- ——- - 
















=R- 








LYNOL g0) 
100 - 064 ¿90€ 4557 





CATA MITRA MA 24 nad AM, 
ahs egress A nr oa DÀ 
a 93 эта м)4105 NO NAL ANO) 1830 
C»5vOQ i75!) 5v 58 NG 87123962 313721.1)2 14 * 
18 m ail 501238802 Sf памі тг 
25906 34 90 FOS 1m INGENOD NO 13954 930105 
Vara ¿Node NO Wainy iy @ 357, 





UWw 0! 1 40) 


— 


envi 


Б0423%407 и3а09 іһемі “0 suda] 


Зәмун 370905 Y 


+ 


E 


13004 ¿rv? > 


ec ace 
03% 
459- 


— 


ва 





о 


yer 
f 
“оэ 





>” of 


rsr | 11 


+81 





50721 405015 


za 


Electronic Logic and Block Diagram of the Honeywell, Inc. 


Signal Conditioner Model Number 552022-1002-100-000. 


Figure B-6 
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Figure B-7 
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Electronic Block of In Situ Measuring System, 
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Figure B-8 Vertical Cross Section of the Current Lines and the Equi- 
potential Lines for the Simple Cell Using the Top and 
Bottom Electrodes as Current Electrodes, 
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Figure B-9 Vertical Cross Section of the Current Lines and the Equi- 
potential Lines for the Modified Cell Using the Outer Two 
Platinum Point Electrodes as Current Electrodes. 
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B. The flow pattern developed in the modified apparatus was similar 
to Figure B-9, 

C. The cell of the redesigned apparatus developed a flow pattern 
similar to Figure 37 in the text. 

D. The flow pattern for the inner corer ring electrode array and the 
microlateralog pad electrode array are presented in the text in Figures 
37 and 38 respectively. The flow pattern for the miniature guard mode 


electrode array is shown in Figure B-10, 


7 


4. Calibration of Equipment. 

A. Using equation 46 of the text and measuring ms resistivity of 
standard solutions with the Beckman Conductivity Meter the cell constants 
for the various electrode configurations of the simple and modified test- 
ing cells were obtained. The various solutions used, with measured resis- 
tivity, resistance and cell constants for two electrode configurations 
are given in Table B-l, Also given are the geometric cell constants 
obtained from equation 46 by dividing the cross sectional area of the cell 
by the potential electrode spacing. 

В. Figure B-11 is the calibration curve for the redesigned testing 
apparatus, This calibration curve was obtained by filling the cell with 
solium chloride solutions of 0.1, 0.15, 0.2, 0.5 and 0,6 normality and 
measuring the signal conditioner output voltage, 

C. Figure 47 in the text shows the calibration curve for the inner 
corer er electrode array, the microlateralog pad electrode array and the 


miniature guard mode electrode array. 
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Figure B-11 Calibration Curve for the Redesigned Apparatus. 
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APPENDIX C 


TEST METHODS AND. DATA 


1. Test Methods Used on Sediments. 

Figure C-1 lists the general properties of the sediments deter- 
mined in this investigation and also identifies the references where the 
procedural method used is outlined. Data obtained by these mom are 
given in Table 3 and Figures 27, 55 and 56 of the text. 

Porosity (n) and void ratio (e) can be determined from the water 
content (w) and specific gravity (G5) of a saturated sediment by the 
following method: | 








and, n = . “е 7 (C2) 
| e+1 ? 
or n = ase 
Е GwWw+l . 
5 
TABLE C-1 


SEDIMENT PROPERTY AND PROCEDURAL METHOD REFERENCE USED IN TESTING. 


Property Procedural Method Reference 

1. Water Content Lambe (46) 

2. Liquid Limit | Olmsted (47) 

3. Plastic Limit Lambe (46) 

4. Specific Gravity Lambe (46) 

5. Grain Size and Hydrometer Lambe (46) 

6. Sediment Classification Lambe (46) 

7. Descriptive Grain Size Measurements Shepard (48) 

| (Boo and 64) s 

8. pH Greweling and Peech (49) 
9. Organic Content Ackroyd (50) 

10. Undrained Shear | : Hansbo (51) 

11. Mineralogy by X-ray Diffraction Cullity (52) 

Technique 


Cl 





2. Experimental Methods and Data Obtained. 
A. Simple Cell Apparatus. 
The porosity was determined in the simple cell by measuring the 
length of the sediment sample in the test cell (L,), the diameter of the 
cell (D.), the weight of soil in the test cell (W.) and the specific grav- 


ity of the soil (6). Then the total volume of the sample (Ул) 15, 
TT 2 
Vip т^, + (C4) 


where: 


Мы is the volume of the soil, 


Y, is the volume of the voids, 


and, | алта Қы . | (С5) 
65 
Then, 
= E C6 
Vs E V. | (C6) 


and the porosity (n) is, 


п = š (C7) 


ү 
Vr 

The procedure described in the text was used to measure the resis- 
tivity of the sediment and interstitial water. Formation Factors were 
Ве by using equation 3 in the text. Table C~2 lists Formation 
Factor - porosity values obtained for various sediments with the simple 
cell apparatus. 

B. Modified Apparatus. 


The porosity was determined in the modified cell by measuring the 


C2 


TABLE С-2 


FORMATION FACTOR - POROSITY DATA OBTAINED WITH THE SIMPLE CELL 
FOR: OTTAWA SAND, GLACIAL SAND (l-a), GLACIAL SAND (1-b) 
GLACIAL SAND (2), NARRAGANSETT BAY SAND, 

PROVIDENCE SILT, KAOLINITE CLAY, AND ILLITE CLAY. 








OTTAWA SAND GLACIAL SAND (1-a) 


C3 


Я Formation | Formation 
: Porosity (4) Fahr Porosity (%) Factor 
40.5 3.54 46.0 3.19 
40.0 3.69 46.7 3.14 
40.4 6.00 46.5 3.23 
40.9 3.74 47.6 3.10 
34.9 4.26 39.5 3.87 
34.0 4.44 39,9 3.83 
35.4 4.44 40.4 3237 
35.4. 4.36 43.9 3.47 
41.1 3.49 39.6 4.04 
41.6 3.50 35.7 4.46 
42.4 3.48 ад 4.71 
до 3.42 29.9 4.80 
34.6 4.24 39.6 4.00 
34.2 4.24 36.7 4.29 
35.3 4.32 33,9 4.46 
36.8 4.20 980.2 4.63 
40,9 3.73 40,5 3,88 
39.7 4.96 
32.3 4.82 
29.9 9.03 GLACIAL SAND (l-b) 
27.3 9.09 Formation 
35.7 4.43 Porosity (%) D tor 
91.2 4.69 
28.7 4.88 46.8 3.14 
26.0 5.06 47.0 3.08 
41.2 3.43 46.9 3.17 
35.9 4.18 46.9 3.11 
33.8 4.90 44.0 Seal 
99%3 4.54 44,6 Jo. 
32.8 4.61 44.6 3.35 
40.7 3,00 43.9 3.18 
34.0 4.26 43.9 3.43 
41.1 3.48 41.4 3.57 
35.0 4.27 39.6 3.80 
41.3 3.96 35.7 4.25 
35.8 4.29 32.3 4.63 
41,9 3.92 31.2 4.80 





TABLE C-2 (CONT.) 


OTTAWA_SAND 


Porosity (4) 


C5 ,» C5 C» C5 C5 C» C5 њ б 
C» E25 ©з єл -J о X t5 єл 
e eee 


DD m о ро м о 


Formation 
Factor 


4.13 
3.55 
4.32 
3.76 
4.05 
4,32 
4.38 
4.44 
3.29 
4.54 


GLACIAL SAND (2) 


Poros ity (%) 


45.6 
41.4 
38.8 
37.6 
36.9 
34.4 
34.0 
45.1 
40.7 
38.2 
36.7 
39.9 
34.3 
43.9 
42.4 
36.4 
35.2 
33.0 
45.1 
41.3 
38.2 
36.6 
35.0 
34.4 
45.5 
40.2 
38.1 
36.6 
35.0 
33.4 


Formation 
Factor 


3.27 
3.90 
3.72 
3.98 
4.20 
4.58 
4.67 
3.47 
3.65 
4.08 
4.29 
4.55 
4.91 
3.62 
3.79 
4.29 
4.50 
4.86 
3.62 
3.54 
4.18 
4.35 
4.48 
4.75 
3.96 
3.93 
4.23 
4.51 
4.74 
9.05 


C4 


GLACIAL SAND (l-b) 


Porosity (%) 


Formation 
Factor 


3.83 
4.16 
4.53 


4.77 


NARRAGANSETT BAY 


STATION 


Porosity (4%) 


40.9 
35.7 
33.5 
32.3 
31.2 
35.7 
34.6 
33.5 
31.2 
39.1 
38.5 
37.0 
34.8 


С 


Formation 
Factor 


4.21 
4.97 
5.50 
9.66 
5.77 
4,94 
9.26 
5.47 
5.73 
4.53 
4.99 
4.73 
5.02 


PROVIDENCE SILT 


Porosity (%) 


Formation 
Factor 


2.58 
2.86 
3.08 
2.89 


KAOLINITE CLAY 


Porosity (%) : 


91,0 
89.7 
88.0 
87.7 
84.0 
78.7 


Formation 
Factor 


1.20 
1.24 
1.25 
1.28 
1.40 
1.53 





GLACIAL SAND (2) ILLITE CLAY 
Formation Formation 

Porosity (%) Factor Porosity (%) Factor 

42.3 3.74 81.0 1.60 

37.7 4.32 80.0 1.63 

2245 5.20 78.7 1.76 

28.7 9.94 

27.3 5.69 

34.6 4.85 

32.3 5.19 

29.9 5.62 

27.3 5,77 

37.9 4.02 i 

/ LUN 
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CS 








height of the sediment column at the end of the test (he) and the weight 
of the soil Ма Then knowing the cross-sectional area of the sample 
(A) and the specific gravity of the soil (С) the total volume (У) can 


be written as: 


Vp = hA , | (C8) 
and, 

у, = Vip - de р | ‚ (C6) 
or, | | | | i 

V =V a “s Е (C9) 

B 1 e 


and dividing by Ур 


ү 
n wag’ = ] . жұт 0 (C10) 
f Y GAh | 
T 5 f 
or, 
ne = 1 - — 2 (C11) 
f 
where, 
С = LR (C12) 
GA ' 


and п; 18 the porosity of the sediment column of height h The porosity 


£° 
at other consolidation pressure height increments can be obtained by using 


equation Cll and the corresponding height recorded at each pressure 


C6 





increment. Then the water content (w) and void ratio (e) can be obtained 


by the following equations: 





N ; (C13) 
: len | 
ата, 
е 
М = с. (C1) 
5 


Electrical resistivity measurements of the sediment were made at 
each pressure increment and Formation Factors'were da by assuming 
the resistivity of the interstitial water was that of the chamber water. 
Table C-3 contains the incremental height of the sediment column and the 
corresponding consolidation pressure, Formation Factor, porosity, void 
ratio, water content, and C-constants (equation C12) for illite and kool- 
inite clays obtained with the modified testing apparatus. 

C. Redesigned Apparatus. 

Porosity was determined at the beginning and at the end of a test 
by measuring the water content and specific gravity of the soil and using 
equation C3. These porosity values were then compared to Formation Factor 
values determined by the method described in the text. Figure 29 in the 
test is an example of the recorded data obtained with the redesigned | 
equipment. If the total test time (ti) is known, along with the time 
prior to rapid decrease in Formation Factor values (t,) and the time For- 
mation Factor values again became relatively constant (t5), and if the 
corresponding Formation Factors at these two times (t, and to) is also 


known, then the recorded curves, such as Figure 29, can be produced, 


C7 


TALLE C-3 


INCREMENTAL SEDIMENT SAMPLE HEIGHT WITH CORRESPONDING CONSOLIDATION 
PRESSURE, FORMATION FACTOR, POROSITY, VOID RATIO, WATER CONTENT, AND 


C-CONSTANT FOR KAOLINITE AND ILLITE CLAY OBTAINED WITH THE 
MODIFIED TESTING APPARATUS. 





(Kaolinite Clay) 








к e 
© с> е 
qe Па ; ES 
4 © ¡er em 
Я 4 TD Ф О м > Ф 
от ін or 44 м2 
B BET oh чо 1 O 
ua ии ү о TD $ 
29 Ж S 3 $E 55 
но O Qs fu fu п. ~ > Са 
20.9 =- 1.99 76.8 3.31 
18.3 .25 l.83 74,0 2.84 
19.9 „40 1„94 70.1 2.34 
14.4 .80 2.13 66.9 2.02 
5,3 1.65 2.42 64.2 1.79 
13.0 2.00 2.54 63.4 1.73 
12.4 3.00 2572 61.6 1.60 
20.3 = 1.56 70.4 3.23 
20.1 „20 1.68 76.1 3.19 
10.9 640 1.81 73.2 2.73 
14,8 . 80 2.05 67.6 2.08 
14.3 2.00 2.14 66.4 15800 
21.2 = 1.64 77.7 3.48 
20.3 529 1.77 76.7 3.29 
18.0 „42 1„83 73.7 2.81 
14.4 „895 2217 67.2 2.04 
13.2 1.65 2.46 64.2 1079 
12.2 3.30 2.79 61.2 1.98 
11.4 6.00 3.02 58.5 1.41 
16.3 = 1„69 74.0 2.85 
14.4 225 1.97 70,6 2.40 
12.9 „40 2.21 67.0 2.03 
11,9 .80 2.95 63.2 4% 72 
10.8 1.65 2.70 60.8 1.55 
10.5 3.00 2.81 59.7 1.48 
16.2 = 1:69 76.5 3.26 
14,1 „25 1.83 73.0 2.71 
12.6 640 2.04 69,8 2.32 
ШІ,1 „80 2.36 65.8 1.92 
10.6 1.00 2.39 64.2 1r 


Content (7) (w) 


Water 


127.18 
109.41 


90.00 
77.88 
69.00 
66,57 
61.72 
124.17 
122.60 
104,95 
80.14 
76.11 
133.93 
126.61 
107.89 
78.62 
68,88 
60.74 
54.24 
109,75 
92.47 
77.92 
66.10 
99.61 
56.92 
125.50 
104.25 
89,07 
73.89 
68.83 





e e ¡C=-Constant 
(cm) (С) 


къ a > 
ч чч 
о оо 





C8 





h Pc FF n e W с 
16.1 с 1.63 74.6 2„94 112.94 4.09 
14.4 „25 1.84 71.6 2:82 96.95 4.09 

12.5 .40 2.14 67.3 2.06 79,09 4.09 
И 25 „80 2.46 63.8 1.76 67.80 4.09 
10,5 1.70 2.68 61.1 1.57 60.28 4.09 
16.8 = 1.67 74.8 2.96 113.93 4.24 
14.7 „25 1.92 1122 2.47 94.88 4.24 
13.2 „40 2208 67.9 2511 81.28 4,24 
12.0 .80 2255 64.7 1.83 70.39 4.24 
11.0 1.65  , 2.84 61.5 1.59 61.32 4.24 
10.5 3.30 2.89 59.6 1.48 56.79 4.24 
15.9 E 1.70 73.4 2.76 106.11 4.23 
18,5 os 1.83 70.8 2.43. 93.38 4.23 
15.2 .40 2.10 68.0 28812 ‘81.68 4923 
11.9 „80 9232 64.5 1.81 69.74 4.23 
11.0 1.65 2.59 61.5 1.60 61.56 4.23 
10.7 2.00 2.65 60.5 1.53 58.83 4.23 
16.2 a 1.65 74.9 2.98 114.62 4.07 
11.8 „80 2.41 65.5 1.90 73.05 4.07 
10.6 1.65 2.79 61.6 1.60 61.71 4.07 
20.6 a 1.80 75.9 378 121.28 4.96 
20.3 #25 1.80 75.6 3.09 118.95 4.96 
15.8 „40 2.03 68.6 2.19 84.05 4.96 
14.4 .80 2.29 65.6 1.90 73.20 4.96 
13.3 1.65 2.56 62.7 1.68 64.67 4.96 
12.9 3.30 2.78 60.3 1252 98.47 4.96 
11.6 6.60 3.10 99.2 1.34 51.49 4.96 

(Illite Clay) 

h Pc FF n е м co 
19.7 ЕЁ 2.05 70.2 2.36 84.59 5.91 
18.8 .40 2.16 68.6 2.18 78.30 5.91 
ay .7 „80 2.28 66,6 1.99 71.47 5.91 
16.9 1.65 2.38 65.0 1.86 66.56 5697 
16.6 2.00 2.50 64.0 1.78 63.72 5.91 
20.1 = 2.12 70.2 2.42 86.60 5.98 
19.5 „40 2.14 69.4 2.29 81.29 5.98 
18.3 ‚85 2,20 67.3 2.06 73.77 5.98 
17.3 1.70 2028 65.4 1.89 67.75 5.98 
16.2 3.30 2.78 63.1 QT 61.29 5.98 
15.9 4.00 2.90 62.4 1.65 59.48 5.98 
19.3 - 2.08 70.8 2.42 86.91 5.64 
18.6 „40 ' 2.20 69.7 2.30 82.45 5.64 
17.1 .70 2.39 67.0 2.03 22.17 5.64 


TABLE C-3 (CONT.) 





(Kaolinite Clay) 








co 





TABLE C-3 (CONT.) 


(Illite Clay) 


h Pc FF n e м 
19.6 - 1.94 70.5 2,39 85.66 
18.8 .20 1299 69,3 2.26 80.91 
18,2 .45 2.20 68.2 2.14 76.87 
17.4 . «09 2.28 66.8 2.01 724513 
16.6 1.65 2.40 65.2 1„87 67.15 
15.6 3.30 2.63 62.9 1.70 60.77 
14.7 6.00 2„93 60.7 1.54 59.36 
18.5 .20 2.09 69.2 2.29 80.69 
16.3 ES 2.16 65.1 1.86 66.83 
15,4 1.65 2,28 63.1 de, 71 61.16 
14.6 3.30 2.44 61.0 1„57 96.13 
14.2 5.90 2.60 99.9 150 93.61 
ЕВ,2 „20 2.14 68.4 2.16 77.41 
ag. 7 .40 2.16 67.4 2.07 74.30 
16.9 . 85 2.24 65.9 1.93 69.32 
16.1 1.65 2.38 64,2 1.80 64.34 
15.2 3.30 2.56 62.1 1.64 58.74 

‚ 14.4 6.00 2.80 60.0 1.50 93.76 
19.5 .20 1.95 69,6 2329 82.22 
19,2 „40 1.98 69.2 2.24 80.40 
17.0 „85 2.05 65.2 1.87 67,08 
16.4 1.65 2.20 63.9 1.74 63.45 
16,2 2.00 2427 63.5 1.74 62.24 
17.7 „20 2.08 67.9 2.14 79.69 
17.3 „40 2„16 67.1 2.04 73.13 
16.5 „85 2421 65.5 1.90 68.10 
15.6 1.65 2.36 63.5 1.74 62.42 
14.7 4.00 2.64 61.3 1.58 56.76 
20.8 229 2.00 68.7 2.15 77.06 

20.7 „40 25127 68.1 2.13 76.34 
20.6 „85 2„14 67.9 24012 79,85 
18.3 1.65 2921 63.9 1777 63.39 
17.2 3.30 2.39 61.6 1.60 57.42 
16.3 6.00 2.63 59.5 1.47 92.09 
20.8 420 2.09 68,4 2.17 77.62 
20.2 „40 2.19 67.5 2.08 74.37 
19.9 .80 2.19 67.0 2.03 72.71 
19.6 1.09 2.40 66.9 1296 71.09 
20.5 „20 1.98 68.2 2.15 77.01 
18.6 .40 2.09 65.0 1.86 66.96 
18.4 1.00 2508 64.6 1.83 65.47 
18.0 2.00 2.26 63.8 1.76 63.25 


Q 
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Table C-4 contains the data obtained with the redesigned equipment, the total 
test time, (t), times ti and tor along with the corresponding Formation 
Factors, porosity, void ratios and water contents for these times. Also 
in this figure is the Yeu of the sediment column at the beginning and at 
the end of a test, the water collected in the reservoir and the consolidation 
pressure used in each test. | 
D, Statistical Technique and Methods. 
Statistical analysis of data was accomplished with a statistical 

ИА of computer programs called **STATPACK (53). The computer programs 
are available at remotely located terminals and prior knowledge of computer 
| ог progranming is not necessary. The analyses of STATPACK eee in this 
study were: 

Transformation to obtain logarithmic functions of Formation 


Factor and porosity data, 


Regression to obtain a functional relationship among For- 
mation Factor (FF) and porosity (n). The rela- 


tionship considered here is of the form: 
Log FF = log a - m log n, 


where FF is called the dependent variable and 
n the independent variable. The slope of the 
line is (m) and the intercept value is (a). 
Line of best fit and the coefficients m and a 
were determined by the method of least squares, 
Other statistics such as, correlation coeffic- 


jent, sum of squares, attributed to the 
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regression, sum of the squares of deviations 
from the regression, F-values for analysis of 
variance, standard error of estimate, standard 
error of regression coefficient, computed t- 


value, and residuals were also obtained, 


Scatter Diagram to obtain Formation Factor — porosity printed 
scatter diagrams of data. The line of best fit 
and the T percent porosity error line were 


/ 


also placed on the scatter diagrams, 


Tables C-5 and C-6 along with Figure C-1 outline the procedure used 
on all sediments analyzed. Table C-5 shows the values of Formation Factor 
porosity data obtained for Providence silt. Also given are the values of 
‘the logarithmic transformation of the Formation Factor and porosity. Table 
С-6 shows the intercept, regression coefficient, standard error of regres- 
sion coefficient, computed t-value, correlation coefficient, standard error 
of estimate and the analyssi of variance for the regression technique. 
Figure C-l is a scatter diagram of the Formation Factor = porosity data. 
The solid line is the line of best fit and the two dashed lines aro the 


Ш 2 percent porosity error lines. 
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Scatter Plot of the Data in Table C-5 with Log Formation 
Factor as the Horizontal Variable and the Log Porosity as 
the Vertical Variable. 
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mitting a value of salinity to be obtained. The bottom 
water salinity obtained in this way was in good agreement 
with the interstitial water salinity of the cores. The 
system detected a core loss, touch down, and sample dis- 
turbance. This relatively low cost system allows accurate, 
quick in situ porosity determination and shows promise in 
emer areas such as the monitoring of pollution, the pre- 
diction of other bottom sediment properties and as a warning 
system on ocean structures. 
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